(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 80122, 1746] NotebookOptionsPosition[ 77166, 1642] NotebookOutlinePosition[ 77572, 1659] CellTagsIndexPosition[ 77529, 1656] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "\:041a\:043e\:044d\:0444\:0444\:0438\:0446\:0438\:0435\:043d\:0442\:044b\ ", " ", "\:043b\:0438\:043d"}], "-", RowBox[{ "\:0439", " ", "\:043c\:043e\:0434\:0435\:043b\:0438", " ", "\:041b\:0410"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"a1", "=", RowBox[{"{", RowBox[{"0.635", ",", "1.38", ",", "2.5"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a2", "=", RowBox[{"{", RowBox[{"5.47", ",", "25.5", ",", "54.4"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a3", "=", RowBox[{"{", RowBox[{"2.72", ",", "9.2", ",", "10.6"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a4", "=", RowBox[{"{", RowBox[{"0.269", ",", "0.527", ",", "0.905"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a5", "=", RowBox[{"{", RowBox[{"3.26", ",", "13.1", ",", "10.4"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a6", "=", RowBox[{"{", RowBox[{"0.709", ",", "1.13", ",", "1.6"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a7", "=", RowBox[{"{", RowBox[{"0.043", ",", "0.0709", ",", "0.0554"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b1", "=", RowBox[{"{", RowBox[{"3.1", ",", "7.32", ",", "12.6"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b2", "=", RowBox[{"{", RowBox[{"20.2", ",", "61.7", ",", "176"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b3", "=", RowBox[{"{", RowBox[{"17.6", ",", "51.2", ",", "33.5"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b4", "=", RowBox[{"{", RowBox[{"0.0719", ",", "0.036", ",", "0.024"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b6", "=", RowBox[{"{", RowBox[{"0.0571", ",", "0.126", ",", "0.178"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b7", "=", RowBox[{"{", RowBox[{"0.0649", ",", "0.0134", ",", "0.0016"}], "}"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.5160849190391197`*^9, 3.5160849588573976`*^9}, { 3.5160850226340456`*^9, 3.516085030441492*^9}, {3.516085068057644*^9, 3.5160851789219847`*^9}, {3.516085217050165*^9, 3.516085973833932*^9}, { 3.516086098680951*^9, 3.516086105950564*^9}, {3.516377782996093*^9, 3.5163777831832933`*^9}, {3.54666247997581*^9, 3.54666248557622*^9}, { 3.5466636866847296`*^9, 3.546663956268803*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ "\:0412\:044b\:0431\:043e\:0440", " ", "\:0440\:0435\:0436\:0438\:043c\:0430", " ", "1"}], ".."}], "3"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"rezh", "=", "2"}], ";"}]}]], "Input", CellChangeTimes->{{3.51608603511084*^9, 3.516086077277714*^9}, { 3.5163747311961317`*^9, 3.5163748030654583`*^9}, {3.5466625023618493`*^9, 3.546662528242295*^9}, {3.5466628391508408`*^9, 3.546662839306841*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0435", " ", "\:0447\:0438\:0441\:043b\:0430", " ", "\:0410\:041f", " ", "\:043a\:0440\:0435\:043d\:0430", " ", "\:0441", " ", "\:0416\:041e\:0421"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tr", "=", "3"}], ";"}], RowBox[{"(*", " ", RowBox[{ "\:0437\:0430\:0434\:0430\:0434\:0438\:043c", " ", "\:0432\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:0438\ \:044f"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"myu", "=", RowBox[{ RowBox[{"(", RowBox[{"18", "-", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "rezh", "]"}], "]"}], "*", "tr"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"b3", "[", RowBox[{"[", "rezh", "]"}], "]"}], "*", "tr"}], ")"}]}]}], "\[IndentingNewLine]", RowBox[{"ii", "=", RowBox[{"108", "/", RowBox[{"(", RowBox[{ RowBox[{"b3", "[", RowBox[{"[", "rezh", "]"}], "]"}], "*", "tr", "*", "tr"}], ")"}]}]}], "\[IndentingNewLine]", RowBox[{"nyu", "=", RowBox[{"216", "/", RowBox[{"(", RowBox[{ RowBox[{"b3", "[", RowBox[{"[", "rezh", "]"}], "]"}], "*", "tr", "*", "tr", "*", "tr"}], ")"}]}]}]}]}]], "Input", CellChangeTimes->{{3.516086085920129*^9, 3.5160860874801316`*^9}, { 3.5160861654178686`*^9, 3.5160861924371157`*^9}, {3.5160862636512413`*^9, 3.5160865646537695`*^9}, 3.5466623708536186`*^9, {3.546662579831586*^9, 3.546662613699245*^9}, {3.5466629220961866`*^9, 3.546662923281789*^9}}], Cell[BoxData[ RowBox[{"-", "0.025781250000000002`"}]], "Output", CellChangeTimes->{ 3.516374492671713*^9, 3.516374629951954*^9, {3.5163747366873417`*^9, 3.5163748067470646`*^9}, 3.516374849007539*^9, 3.516375498995881*^9, 3.516375538838351*^9, 3.5163758644577227`*^9, {3.5163759479802694`*^9, 3.516375966622302*^9}, 3.5163759995851603`*^9, 3.5163765214060764`*^9, 3.5163765928854017`*^9, 3.5163768389666348`*^9, 3.5163769426372166`*^9, 3.516377036174981*^9, 3.516377259192973*^9, {3.516377330063897*^9, 3.516377352590337*^9}, 3.516377493318184*^9, 3.5163775272170434`*^9, 3.5163775875423493`*^9, 3.5163776216596093`*^9, {3.516377788814903*^9, 3.516377803213728*^9}, 3.5163779144731236`*^9, {3.516378053843768*^9, 3.5163780836866207`*^9}, 3.516417858520227*^9, {3.516417949125186*^9, 3.5164179649904137`*^9}, 3.5164180473273582`*^9, 3.516591173138672*^9, 3.546662413114093*^9, {3.546662512533067*^9, 3.546662533405904*^9}, { 3.5466625905644045`*^9, 3.546662618363653*^9}, 3.54666284451725*^9, 3.5466628882597275`*^9, 3.546662926214594*^9, 3.54666396018441*^9, 3.546664110989875*^9, 3.546665023061077*^9, 3.5466653080735774`*^9, 3.546665370582887*^9, 3.546666687423608*^9, 3.54667156222937*^9, 3.5466723265527124`*^9, 3.5467439322050533`*^9, 3.5467439684439173`*^9, 3.546924846115532*^9, 3.546924912805649*^9, 3.546925008449417*^9}], Cell[BoxData["0.23437499999999997`"], "Output", CellChangeTimes->{ 3.516374492671713*^9, 3.516374629951954*^9, {3.5163747366873417`*^9, 3.5163748067470646`*^9}, 3.516374849007539*^9, 3.516375498995881*^9, 3.516375538838351*^9, 3.5163758644577227`*^9, {3.5163759479802694`*^9, 3.516375966622302*^9}, 3.5163759995851603`*^9, 3.5163765214060764`*^9, 3.5163765928854017`*^9, 3.5163768389666348`*^9, 3.5163769426372166`*^9, 3.516377036174981*^9, 3.516377259192973*^9, {3.516377330063897*^9, 3.516377352590337*^9}, 3.516377493318184*^9, 3.5163775272170434`*^9, 3.5163775875423493`*^9, 3.5163776216596093`*^9, {3.516377788814903*^9, 3.516377803213728*^9}, 3.5163779144731236`*^9, {3.516378053843768*^9, 3.5163780836866207`*^9}, 3.516417858520227*^9, {3.516417949125186*^9, 3.5164179649904137`*^9}, 3.5164180473273582`*^9, 3.516591173138672*^9, 3.546662413114093*^9, {3.546662512533067*^9, 3.546662533405904*^9}, { 3.5466625905644045`*^9, 3.546662618363653*^9}, 3.54666284451725*^9, 3.5466628882597275`*^9, 3.546662926214594*^9, 3.54666396018441*^9, 3.546664110989875*^9, 3.546665023061077*^9, 3.5466653080735774`*^9, 3.546665370582887*^9, 3.546666687423608*^9, 3.54667156222937*^9, 3.5466723265527124`*^9, 3.5467439322050533`*^9, 3.5467439684439173`*^9, 3.546924846115532*^9, 3.546924912805649*^9, 3.546925008449417*^9}], Cell[BoxData["0.15625`"], "Output", CellChangeTimes->{ 3.516374492671713*^9, 3.516374629951954*^9, {3.5163747366873417`*^9, 3.5163748067470646`*^9}, 3.516374849007539*^9, 3.516375498995881*^9, 3.516375538838351*^9, 3.5163758644577227`*^9, {3.5163759479802694`*^9, 3.516375966622302*^9}, 3.5163759995851603`*^9, 3.5163765214060764`*^9, 3.5163765928854017`*^9, 3.5163768389666348`*^9, 3.5163769426372166`*^9, 3.516377036174981*^9, 3.516377259192973*^9, {3.516377330063897*^9, 3.516377352590337*^9}, 3.516377493318184*^9, 3.5163775272170434`*^9, 3.5163775875423493`*^9, 3.5163776216596093`*^9, {3.516377788814903*^9, 3.516377803213728*^9}, 3.5163779144731236`*^9, {3.516378053843768*^9, 3.5163780836866207`*^9}, 3.516417858520227*^9, {3.516417949125186*^9, 3.5164179649904137`*^9}, 3.5164180473273582`*^9, 3.516591173138672*^9, 3.546662413114093*^9, {3.546662512533067*^9, 3.546662533405904*^9}, { 3.5466625905644045`*^9, 3.546662618363653*^9}, 3.54666284451725*^9, 3.5466628882597275`*^9, 3.546662926214594*^9, 3.54666396018441*^9, 3.546664110989875*^9, 3.546665023061077*^9, 3.5466653080735774`*^9, 3.546665370582887*^9, 3.546666687423608*^9, 3.54667156222937*^9, 3.5466723265527124`*^9, 3.5467439322050533`*^9, 3.5467439684439173`*^9, 3.546924846115532*^9, 3.546924912805649*^9, 3.546925008449417*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041c\:0430\:0442\:0440\:043c\:0446\:0430", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0441\:0432\:043e\:0431\:043e\:0434\:043d\:043e\:0433\:043e", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"MurKrap0", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", RowBox[{"b1", "[", RowBox[{"[", "rezh", "]"}], "]"}]}], ")"}], ",", "0", ",", RowBox[{"b3", "[", RowBox[{"[", "rezh", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myu", ",", RowBox[{"nyu", "/", "s"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0441", " ", "\:0437\:0430\:0434\:0430\:043d\:043d\:044b\:043c", " ", "\:0443\:0433\:043b\:043e\:043c", " ", "\:043a\:0440\:0435\:043d\:0430"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"MurKrap", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", RowBox[{"b1", "[", RowBox[{"[", "rezh", "]"}], "]"}]}], ")"}], ",", "0", ",", RowBox[{"b3", "[", RowBox[{"[", "rezh", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "s", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myu", ",", RowBox[{"(", RowBox[{"ii", "+", RowBox[{"nyu", "/", "s"}]}], ")"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.516086623450273*^9, 3.516086650875121*^9}, 3.5160876413072233`*^9, {3.5160876946437173`*^9, 3.5160877272009745`*^9}, { 3.5160884451454353`*^9, 3.516088454708252*^9}, {3.51608849320912*^9, 3.516088511663952*^9}, 3.5160893877642913`*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"7.32`", "\[VeryThinSpace]", "+", "s"}], "0", "51.2`"}, { RowBox[{"-", "1"}], "0", "0"}, { RowBox[{"-", "0.025781250000000002`"}], FractionBox["0.15625`", "s"], RowBox[{"-", "1"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{ 3.5163744927965136`*^9, 3.5163746300143547`*^9, {3.516374736796542*^9, 3.516374806871865*^9}, 3.516374849069939*^9, 3.5163754990426807`*^9, 3.516375538900751*^9, 3.516375864504523*^9, {3.516375948027069*^9, 3.5163759666691017`*^9}, 3.51637599963196*^9, 3.5163765214528766`*^9, 3.5163765929322023`*^9, 3.5163768390290346`*^9, 3.5163769426840167`*^9, 3.516377036221781*^9, 3.516377259239773*^9, {3.5163773301106973`*^9, 3.5163773526371365`*^9}, 3.516377493364984*^9, 3.5163775272638435`*^9, 3.5163775875891495`*^9, 3.5163776217064095`*^9, {3.516377788939703*^9, 3.5163778032605286`*^9}, 3.5163779145199237`*^9, {3.5163780539061685`*^9, 3.516378083749021*^9}, 3.516417858723027*^9, {3.516417949171986*^9, 3.5164179650372143`*^9}, 3.516418047374159*^9, 3.516591173357422*^9, 3.5466624132388926`*^9, {3.546662512564267*^9, 3.546662533437104*^9}, { 3.5466625905956044`*^9, 3.5466626183792534`*^9}, 3.5466628445328503`*^9, 3.546662888290927*^9, 3.5466629262457943`*^9, 3.5466639602000103`*^9, 3.546664111021075*^9, 3.546665023076677*^9, 3.5466653081047773`*^9, 3.546665370598487*^9, 3.5466666874548073`*^9, 3.54667156224497*^9, 3.5466723265839124`*^9, 3.546743932251854*^9, 3.5467439684595175`*^9, 3.546924846333932*^9, 3.5469249128680487`*^9, 3.546925008496217*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"7.32`", "\[VeryThinSpace]", "+", "s"}], "0", "51.2`"}, { RowBox[{"-", "1"}], "s", "0"}, { RowBox[{"-", "0.025781250000000002`"}], RowBox[{"0.23437499999999997`", "\[VeryThinSpace]", "+", FractionBox["0.15625`", "s"]}], RowBox[{"-", "1"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{ 3.5163744927965136`*^9, 3.5163746300143547`*^9, {3.516374736796542*^9, 3.516374806871865*^9}, 3.516374849069939*^9, 3.5163754990426807`*^9, 3.516375538900751*^9, 3.516375864504523*^9, {3.516375948027069*^9, 3.5163759666691017`*^9}, 3.51637599963196*^9, 3.5163765214528766`*^9, 3.5163765929322023`*^9, 3.5163768390290346`*^9, 3.5163769426840167`*^9, 3.516377036221781*^9, 3.516377259239773*^9, {3.5163773301106973`*^9, 3.5163773526371365`*^9}, 3.516377493364984*^9, 3.5163775272638435`*^9, 3.5163775875891495`*^9, 3.5163776217064095`*^9, {3.516377788939703*^9, 3.5163778032605286`*^9}, 3.5163779145199237`*^9, {3.5163780539061685`*^9, 3.516378083749021*^9}, 3.516417858723027*^9, {3.516417949171986*^9, 3.5164179650372143`*^9}, 3.516418047374159*^9, 3.516591173357422*^9, 3.5466624132388926`*^9, {3.546662512564267*^9, 3.546662533437104*^9}, { 3.5466625905956044`*^9, 3.5466626183792534`*^9}, 3.5466628445328503`*^9, 3.546662888290927*^9, 3.5466629262457943`*^9, 3.5466639602000103`*^9, 3.546664111021075*^9, 3.546665023076677*^9, 3.5466653081047773`*^9, 3.546665370598487*^9, 3.5466666874548073`*^9, 3.54667156224497*^9, 3.5466723265839124`*^9, 3.546743932251854*^9, 3.5467439684595175`*^9, 3.546924846333932*^9, 3.5469249128680487`*^9, 3.546925008511817*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:0420\:0430\:0441\:0447\:0435\:0442", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:043e\:0439\ ", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:0438"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"w", "=", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"Det", "[", "MurKrap0", "]"}], "/", RowBox[{"Det", "[", "MurKrap", "]"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.51608775893143*^9, 3.5160877890550833`*^9}, 3.516087819646737*^9, {3.516089305552147*^9, 3.5160893120729585`*^9}, { 3.5163733133390303`*^9, 3.5163733494530935`*^9}, 3.516373385661757*^9}], Cell[BoxData[ RowBox[{"-", FractionBox["8.`", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "2.`"}], "-", RowBox[{"1.`", " ", "s"}]}], ")"}], "3"]]}]], "Output", CellChangeTimes->{ 3.5163744928589134`*^9, 3.5163746300611544`*^9, {3.516374736858942*^9, 3.5163748069498653`*^9}, 3.5163748491167393`*^9, 3.516375499120681*^9, 3.516375538947551*^9, 3.516375864551323*^9, {3.5163759480738697`*^9, 3.516375966700302*^9}, 3.51637599967876*^9, 3.5163765214996767`*^9, 3.516376592979002*^9, 3.5163768390758348`*^9, 3.516376942730817*^9, 3.5163770362685814`*^9, 3.5163772592865725`*^9, {3.5163773301574974`*^9, 3.516377352683937*^9}, 3.516377493411784*^9, 3.5163775273106437`*^9, 3.5163775876359496`*^9, 3.5163776217532096`*^9, {3.516377789017703*^9, 3.516377803307328*^9}, 3.516377914566724*^9, {3.5163780539373684`*^9, 3.516378083795821*^9}, 3.5164178593938284`*^9, {3.5164179492187862`*^9, 3.516417965084014*^9}, 3.5164180474209585`*^9, 3.516591173669922*^9, 3.5466624136132936`*^9, {3.5466625125798674`*^9, 3.5466625334683037`*^9}, { 3.546662590611204*^9, 3.546662618410453*^9}, 3.5466628445640507`*^9, 3.546662888306527*^9, 3.5466629262769938`*^9, 3.5466639602312098`*^9, 3.5466641110366745`*^9, 3.546665023107877*^9, 3.5466653081203775`*^9, 3.5466653706296873`*^9, 3.5466666874704075`*^9, 3.54667156227617*^9, 3.5466723266151123`*^9, 3.5467439322986536`*^9, 3.546743968490717*^9, 3.546924846926733*^9, 3.546924912914849*^9, 3.5469250085586166`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Wt", "=", RowBox[{ RowBox[{"Det", "[", "MurKrap0", "]"}], "/", RowBox[{"Det", "[", "MurKrap", "]"}]}]}]], "Input", CellChangeTimes->{{3.516374249442481*^9, 3.51637427453653*^9}, 3.516374607909116*^9}], Cell[BoxData[ RowBox[{"-", FractionBox["8.`", RowBox[{"s", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "12.`"}], "-", FractionBox["8.`", "s"], "-", RowBox[{"6.`", " ", "s"}], "-", SuperscriptBox["s", "2"]}], ")"}]}]]}]], "Output", CellChangeTimes->{ 3.5163744929057136`*^9, 3.5163746301079545`*^9, {3.5163747369057417`*^9, 3.5163748069810653`*^9}, 3.516374849147939*^9, 3.516375499167481*^9, 3.5163755389787507`*^9, 3.516375864582523*^9, {3.5163759481050696`*^9, 3.5163759667471027`*^9}, 3.51637599970996*^9, 3.516376521530877*^9, 3.5163765930102024`*^9, 3.5163768391070347`*^9, 3.516376942762017*^9, 3.516377036299781*^9, 3.516377259317773*^9, {3.5163773301886973`*^9, 3.5163773527307367`*^9}, 3.516377493458584*^9, 3.5163775273418436`*^9, 3.5163775876671495`*^9, 3.5163776217844095`*^9, {3.5163777890645037`*^9, 3.5163778033385286`*^9}, 3.516377914597924*^9, {3.516378053968569*^9, 3.516378083827021*^9}, 3.516417859503029*^9, {3.5164179492655864`*^9, 3.5164179651152143`*^9}, 3.516418047467759*^9, 3.516591173716797*^9, 3.5466624136444936`*^9, {3.5466625126110673`*^9, 3.546662533499504*^9}, { 3.5466625906268044`*^9, 3.546662618426053*^9}, 3.5466628445796504`*^9, 3.5466628883221273`*^9, 3.546662926292594*^9, 3.54666396024681*^9, 3.5466641110522747`*^9, 3.546665023123477*^9, 3.546665308151578*^9, 3.5466653706452875`*^9, 3.546666687501608*^9, 3.54667156229177*^9, 3.5466723266307125`*^9, 3.546743932314254*^9, 3.546743968506317*^9, 3.546924846942333*^9, 3.5469249129460487`*^9, 3.546925008589817*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:0435\:0440\:0435\:0439\:0434\:0435\:043c", " ", "\:043e\:0442", " ", "s", " ", "\:043a", " ", "\:043a\:043e\:043c\:043f\:043b\:0435\:043a\:0441\:043d\:043e\:0439", " ", "\:0444\:043e\:0440\:043c\:0435"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"wk", "=", RowBox[{"w", "/.", RowBox[{"s", "\[Rule]", RowBox[{"I", "*", "omega"}]}]}]}]}]], "Input", CellChangeTimes->{{3.5160879222637167`*^9, 3.5160879667549953`*^9}, { 3.516088011917074*^9, 3.5160880449579325`*^9}}], Cell[BoxData[ RowBox[{"-", FractionBox["8.`", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "2.`"}], "-", RowBox[{ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"1.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "omega"}]}], ")"}], "3"]]}]], "Output", CellChangeTimes->{ 3.5163744929369135`*^9, 3.5163746301547546`*^9, {3.516374736936942*^9, 3.5163748070278654`*^9}, 3.5163748491947393`*^9, 3.516375499214281*^9, 3.516375539025551*^9, 3.516375864629323*^9, {3.51637594815187*^9, 3.516375966778302*^9}, 3.5163759997567606`*^9, 3.5163765215620766`*^9, 3.516376593057002*^9, 3.516376839138235*^9, 3.5163769427932167`*^9, 3.5163770363309813`*^9, 3.516377259364573*^9, {3.5163773302198973`*^9, 3.516377352761937*^9}, 3.5163774934897842`*^9, 3.5163775273886433`*^9, 3.5163775877139497`*^9, 3.5163776218312097`*^9, {3.5163777891113033`*^9, 3.5163778033853283`*^9}, 3.516377914644724*^9, {3.5163780540153685`*^9, 3.5163780838738213`*^9}, 3.516417859534229*^9, {3.5164179492967863`*^9, 3.5164179651464148`*^9}, 3.516418047514559*^9, 3.516591173763672*^9, 3.5466624136600933`*^9, {3.5466625126266675`*^9, 3.5466625335307045`*^9}, { 3.5466625906580043`*^9, 3.5466626184416533`*^9}, 3.5466628446108503`*^9, 3.5466628883533278`*^9, 3.546662926308194*^9, 3.54666396026241*^9, 3.546664111083475*^9, 3.5466650231390767`*^9, 3.5466653081671777`*^9, 3.5466653706608872`*^9, 3.546666687517208*^9, 3.54667156230737*^9, 3.5466723266619124`*^9, 3.5467439323454537`*^9, 3.5467439685219173`*^9, 3.546924846973533*^9, 3.546924912961649*^9, 3.546925008621017*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"f", "=", RowBox[{ RowBox[{"(", RowBox[{"Abs", "[", "wk", "]"}], ")"}], "^", "2"}]}]], "Input", CellChangeTimes->{{3.5160881542985244`*^9, 3.5160881791025677`*^9}, { 3.5160882708931293`*^9, 3.5160882833419514`*^9}}], Cell[BoxData[ FractionBox["64.`", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"-", "2.`"}], "-", RowBox[{ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"1.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "omega"}]}], "]"}], "6"]]], "Output", CellChangeTimes->{ 3.5163744929837136`*^9, 3.5163746302015543`*^9, {3.516374736983742*^9, 3.516374807074665*^9}, 3.5163748492415395`*^9, 3.516375499261081*^9, 3.5163755390723515`*^9, 3.516375864660523*^9, {3.5163759481830697`*^9, 3.5163759668251023`*^9}, 3.51637599978796*^9, 3.5163765216088767`*^9, 3.5163765930882025`*^9, 3.5163768391850348`*^9, 3.5163769428400173`*^9, 3.516377036377781*^9, 3.516377259395773*^9, {3.516377330266698*^9, 3.5163773528087373`*^9}, 3.516377493520984*^9, 3.5163775274198437`*^9, 3.5163775877451496`*^9, 3.5163776218624096`*^9, {3.516377789158104*^9, 3.5163778034165287`*^9}, 3.5163779146759243`*^9, {3.516378054046569*^9, 3.516378083905021*^9}, 3.516417859612229*^9, {3.5164179493279867`*^9, 3.5164179651932144`*^9}, 3.516418047561359*^9, 3.516591173810547*^9, 3.5466624136912937`*^9, {3.5466625126422677`*^9, 3.546662533561904*^9}, { 3.5466625906736045`*^9, 3.546662618472853*^9}, 3.5466628446264505`*^9, 3.5466628883689275`*^9, 3.546662926339394*^9, 3.54666396029361*^9, 3.5466641110990753`*^9, 3.546665023154677*^9, 3.5466653081827774`*^9, 3.546665370676487*^9, 3.546666687532808*^9, 3.54667156232297*^9, 3.546672326677512*^9, 3.5467439323610535`*^9, 3.5467439685531178`*^9, 3.546924846989133*^9, 3.5469249129928493`*^9, 3.546925008652217*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{"f", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "5"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.51608832249802*^9, 3.516088395911749*^9}, 3.516088542208806*^9, 3.5160886380709743`*^9, {3.5160888670325766`*^9, 3.5160889096206512`*^9}, {3.516088962567144*^9, 3.5160889928155975`*^9}, { 3.516089595806257*^9, 3.516089623199905*^9}, {3.51634370432053*^9, 3.516343706348533*^9}, {3.546743927774646*^9, 3.5467439287418475`*^9}, { 3.5469254976194763`*^9, 3.546925511737501*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0nk81PkfB3BX2ZDw3UNNoX6Vkmy22l2WXp9NbRH9Otj81B6GORAhR6Eo WWwTstJsw7qi5JYUK+vscOWatDuj5Mpi5usa9+Y3+348Xo/X4/nf64/3euaZ YywlBQUFnjz/dvJF5nJT1+A9j5b+PRp1HQa2Jwz+h8YEBbtFuZeeVIYbGLjB et83u8bl1owvLVR75QvqqOjBa7mTq67pLzM4j8bO1rY6uY1p15glTigYVXui MuTu1vzt/WLXZdj+ZfVnkNwbs1Z8Pr8Qgf9UL6u3kfu0pb/njH40rgwVHqLk vt/Zc3vKioeDt855dr6nMe9hKx7nxCJVlmgUI3eqojh6Zv91lEptw4nc+5qv b5ruiofH4d7AkX9oDPEPVE+5JeB8UeYHsXLzXP85NblwA2pRo7u3yd2xwE0Y 0/8VP3qGcOwWaQQ81dtBF90CJzD+dMcCjTW/dDZKrJIg/bt+1TG5mUZEeYST gpXjbp+ZzdNQWgh6MRWSCsX33cVpczQ6exzWTe1PgwUzoENJ7nO56mUTXekI +5s1kz9DwyZ+cPlEagaaeqcX5qZprD1XbT/udhsnu/UzLOWusgocoxcyIdSd FeVP0VAV9RpK9LMR7e2tYTZO41XVY//RoWwwK4m59RiNe1n82pGie+h1FGod o2kc9j38w7BVLl5Xvv3dTkIjcUV54jtOAcKvMTZKhmiUBNr4ntUtxPpjZ4tr 3tFIGv5i9WRIIQTu7uNxgzS8WnU4E/uLkPVZRsQn/TR0kp4pjXcVY/vHuYbN r2ksaJZme1vch3uOxUluN43+SxlHxlLvQ/rLxMFFEY1SzsXfaLcSbDVe5aj5 Jw2nnbvMpQsPwMr62kPUTiOtIcVnRL8MwaFluWvqaUzxHxi94ZXBZV8MY00t jW84jX3tc2VYa3NhjXY1jVHlGYfyjnKkap9tG6ig8aXlYfPoyAqYLBi5GJXQ aC1YVDKk/8Bn/13cviFNvhcrvl9pUgVvY9NHgt9ofLU5wTzStQrrfhXMaSTR iJnKngxqq4Lf0eV5LxNp7LouZDFzqpGdWbNRlUcjtMHYxvSHWmTqDL577Ufj btHDzTk3avFWkgKhD412/l7lTU21EK01vVjrJf9fjmOFrlkd2Hd7KyK5NBqU r5goUPUYvPhTcsFJGh9binRePHmCGkeGNJTQyC+IFp02aUD3lruBFio0hnWG eTbODTBJ0tIIVqBhGGCzZ0tCA9zeCGaLF6VItVBP65trgLVjyztNmRTxz3ls p/pGUAMPQwMHpAjoix0/cKoZzPybBZr1UqSoN4Zr729FftHWUa9QKfZYBp8K +74V+mGqP78JkkLstW03HdiKO53cGesAKXTbrw423WvFF86DRzQ8pbjOt7WO 0mrDp2cYPgecpAjf1KK5JGrDRTsnVcFOKTho+3XUtwPmEw+9hW8lUEifLhdd 7QDl8FyYLpaAr7JW3HC7A0GPqSb3LgmePWfrZb/swIhv6IXhJgm22C+ms8w7 8ejz5NPZDyUY4hrmvlYU4oZ9vd1WngTu8RcqW+NeQsCjrI+YSuA1sLWvJP9P jCH/mKXnKNTWXzF882UPdIfOq2zgDSNGq9Eo7cc+SFVsbzzkD4Epsv9DuWkQ XJ/tpdWnBuFTf4V5xnAYDNO1YcwN/Sjy1jiuay2B327XHLNPepEW7DG848QY cq8ej31Ev8HllMWjU6PjmAphMGa/FEOho4dUBE2iLeCOvqL4FTb7sc/ZG8kw +1FTpuyqEP2ym9+qP55G/48R2dpP22Ez4H4h3mcWO0P8lvnoteJwf9w2Pb15 cI3rhLzcRgzXMsOe5S1AzWdDkq/sKUy0X66z+uEfjEiOX6oIqIMaZM7iyfeQ GmWmRGdV44DEaWaDgwLZNX1a0KpViVc1ropOxoqkyv6rjriCMjyIqvveQKZI tsyP3VixoxTsZf49LS1K5Lg738xlohguvsEPNG4pk+qd55+1fFeITznbWmyd VUhY8r2Mn11y8ULHWN3fbBnplB075HnqLn7PvywwVl5OdnVVhRZvzkSwotGR pu7lxEd17OBKk3SoiYMOepWoko96a4V5rBQMZK1PSAn5gOReO+Bf6ClAwpbV KmUOK4jv5B+P7StuIu4W+6mzgRqJcVYVTz9JgFPCdy0a79SIie2Qzfl98Zjd lzmXUqlOhDd77zcJY9E3sPdcTaQGAfOKZ7YGD+mWpoxm55UkvT8f8WFRyIvw 1+Zs1yTfOJ7Q2pQRgUgHl0PKMk1SkcSvZ+wIx25xWnZk8yoSzr8V6dkTBv5W QV7zIS0y+BW7+cTSBYRWXeMfKtMimyRaLz44GwQ7hkNmw0Zt8ii3fP01j0B8 7fX0pHa0NtFSF5wUjPshZS1zY5NMm5T4fHvDL9EXG2vvJB4+qUNWx1U+kt3x Rtz2GqP2Sh1yM3TETX/aE+xJlbbnRhTZnGrk5FPlDmbRhEF0LEV0BQ6X5tdw cfry8UWb6xQp9zNuF63iwv94SZfGLxTxzVw9X6PCRZTMPyYukSKqEXHV6VIO 8s3mFm4mU0TNW/ziTg0Hc9VLL7NyKKJ5wGW2zoODuHaNa3VPKXI5ye7Urho2 bmV4cn96TpGfRwsZDQ/ZyPBrsTrYSJGxESVrVh4bpR/HzTe2UGTgyKueHD4b YieK2yGkSOT8hbzzZ9gw7Ftt1dtHkd/LXyt9pMfGjpIgvdsDFHE027pN4UM2 zCNEc6x3FHE17I4dW8GGrWFy4d/DFNGycFV6K2PB18NAb3ycIsrdMV8otLAQ YnFprniSIn+NlEo/rGMhYmVvp5+MIqN/ocGknAV+QcbV2VmKqGsn9PpksZAe psIpn6dIU//Z1YIkFnKOsvaGLFJEUNDNfR7PQsmGJ+v2vKcIw6WneTGKhcrJ zXNLSxSpa+Lt3x3Kwv8BBL1s2w== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, ImageSize->{546., Automatic}, PlotRange->{{0, 5}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Automatic}]], "Output", CellChangeTimes->{ 3.516374493108514*^9, 3.516374630326355*^9, {3.5163747371241426`*^9, 3.5163748071994658`*^9}, 3.51637484938194*^9, 3.516375499417082*^9, 3.516375539197151*^9, 3.516375864785323*^9, {3.51637594830787*^9, 3.5163759669343023`*^9}, 3.5163759999439607`*^9, 3.516376521733677*^9, 3.516376593197403*^9, 3.5163768393098354`*^9, 3.516376942949217*^9, 3.5163770364869814`*^9, 3.5163772595049734`*^9, {3.5163773303914976`*^9, 3.5163773529179373`*^9}, 3.5163774936457844`*^9, 3.5163775275446444`*^9, 3.5163775878699503`*^9, 3.5163776219872103`*^9, {3.516377789298504*^9, 3.516377803541329*^9}, 3.5163779147851243`*^9, {3.5163780541557693`*^9, 3.516378084014221*^9}, 3.51641786009583*^9, {3.5164179494527864`*^9, 3.516417965302415*^9}, 3.516418047686159*^9, 3.516591174169922*^9, 3.5466624141904945`*^9, {3.5466625127202673`*^9, 3.546662533655504*^9}, { 3.546662590751605*^9, 3.5466626185352535`*^9}, 3.5466628447044506`*^9, 3.5466628884625273`*^9, 3.546662926417394*^9, 3.5466639603560104`*^9, 3.546664111161475*^9, 3.5466650232170773`*^9, 3.5466653082451773`*^9, 3.5466653707544875`*^9, 3.546666687595208*^9, 3.5466715623853703`*^9, 3.5466723267555127`*^9, 3.546743932563854*^9, 3.5467439685687175`*^9, 3.546924847956335*^9, 3.546924913039649*^9, 3.5469250086990175`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"p", "=", RowBox[{"NIntegrate", "[", RowBox[{"f", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5160898890711718`*^9, 3.5160899559640894`*^9}, { 3.516343711824143*^9, 3.516343713181345*^9}, {3.546924895754819*^9, 3.5469248965192204`*^9}}], Cell[BoxData["1.1780972450963092`"], "Output", CellChangeTimes->{3.5469249133672495`*^9, 3.5469250087614174`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NIntegrate", "[", RowBox[{"f", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "\[Infinity]"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5469249007936277`*^9, 3.546924901495629*^9}}], Cell[BoxData["1.1780972450963092`"], "Output", CellChangeTimes->{3.54692491341405*^9, 3.546925008808217*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"omega0", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"shag", "=", "0.1"}], ";"}], "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{"shag", "\[GreaterEqual]", RowBox[{"10", "^", RowBox[{"-", "5"}]}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"NIntegrate", "[", RowBox[{"f", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "omega0"}], "}"}]}], "]"}], "-", RowBox[{"0.95", "*", "p"}]}], ")"}], "<", "0"}], ",", RowBox[{"omega0", "=", RowBox[{"omega0", "+", "shag"}]}]}], "]"}], ",", RowBox[{"omega0", "=", RowBox[{"omega0", "-", "shag"}]}], ",", RowBox[{"shag", "=", RowBox[{"shag", "/", "10"}]}]}], "}"}]}], "]"}], "\[IndentingNewLine]", "omega0"}], "Input", CellChangeTimes->{{3.5163435178690023`*^9, 3.5163436794072857`*^9}, { 3.5163437220577607`*^9, 3.516343761947031*^9}, {3.5163438643456106`*^9, 3.5163438797740383`*^9}, {3.516372786427702*^9, 3.516372865644641*^9}, { 3.546664995885829*^9, 3.5466650040914435`*^9}, {3.546667828659212*^9, 3.546667862464472*^9}, {3.5469249451913056`*^9, 3.546924945986907*^9}}], Cell[BoxData["2.2991899999999985`"], "Output", CellChangeTimes->{{3.546667851357252*^9, 3.546667863946474*^9}, 3.5466715629937716`*^9, 3.5466723273639135`*^9, 3.5467439334374557`*^9, 3.546743969161518*^9, 3.546924848736336*^9, 3.54692491358565*^9, 3.5469250099782195`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"p95", "=", RowBox[{"NIntegrate", "[", RowBox[{"f", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "omega0"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5163436706244707`*^9, 3.5163436719036727`*^9}, { 3.5163728835690727`*^9, 3.5163729252991457`*^9}, 3.516373024967721*^9, { 3.5469249386392937`*^9, 3.5469249393412952`*^9}}], Cell[BoxData["1.119191720270109`"], "Output", CellChangeTimes->{ 3.516374493826115*^9, 3.516374631028356*^9, {3.5163747378261433`*^9, 3.516374807901467*^9}, 3.516374850068341*^9, 3.5163755001346827`*^9, 3.5163755399147525`*^9, 3.5163758654717245`*^9, {3.516375948978671*^9, 3.5163759676051035`*^9}, 3.516376000661562*^9, 3.516376522404478*^9, 3.516376593883804*^9, 3.5163768399962363`*^9, 3.516376943635618*^9, 3.5163770371733828`*^9, 3.5163772601913743`*^9, {3.516377331062299*^9, 3.5163773536043386`*^9}, 3.5163774943009853`*^9, 3.5163775282154455`*^9, 3.5163775885407515`*^9, 3.516377622658011*^9, {3.516377789953705*^9, 3.51637780421213*^9}, 3.5163779154559255`*^9, {3.5163780548265705`*^9, 3.5163780846694226`*^9}, 3.516417861328232*^9, {3.5164179500143876`*^9, 3.5164179658952155`*^9}, 3.5164180482945604`*^9, 3.516591175123047*^9, 3.5466624175289*^9, {3.5466625132974687`*^9, 3.5466625343887053`*^9}, { 3.546662591406806*^9, 3.5466626190344543`*^9}, 3.5466628452036514`*^9, 3.5466628889773283`*^9, 3.546662927041395*^9, 3.546663960980011*^9, 3.546664111785476*^9, 3.5466650238410783`*^9, 3.546665308869179*^9, 3.5466653713784885`*^9, 3.546666688203609*^9, 3.546671563024971*^9, 3.546672327395114*^9, 3.5467439335154557`*^9, 3.5467439691927185`*^9, 3.546924848767536*^9, 3.54692491361685*^9, 3.5469250100094194`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:0415\:0420\:0418\:041e\:0414", " ", "\:041a\:0412\:0410\:041d\:0422\:041e\:0412\:0410\:041d\:0418\:042f"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Tkvant", "=", RowBox[{"Pi", "/", "omega0"}]}]}]], "Input", CellChangeTimes->{{3.516373076057811*^9, 3.5163731523195453`*^9}, { 3.516591167576172*^9, 3.516591168841797*^9}, {3.5466665777398148`*^9, 3.5466665779426155`*^9}, {3.5466666748811855`*^9, 3.546666680044795*^9}}], Cell[BoxData["1.3663910566720432`"], "Output", CellChangeTimes->{ 3.516374493857315*^9, 3.516374631075156*^9, {3.516374737872944*^9, 3.516374807948267*^9}, 3.516374850115141*^9, 3.516375500165883*^9, 3.5163755399615526`*^9, 3.5163758655029244`*^9, {3.516375949025471*^9, 3.516375967651904*^9}, 3.516376000708362*^9, 3.5163765224512787`*^9, 3.516376593930604*^9, 3.5163768400274363`*^9, 3.5163769436824183`*^9, 3.516377037220183*^9, 3.516377260222574*^9, {3.516377331109099*^9, 3.5163773536511383`*^9}, 3.516377494347786*^9, 3.516377528262245*^9, 3.516377588571951*^9, 3.516377622704811*^9, {3.516377790000505*^9, 3.51637780424333*^9}, 3.5163779155027256`*^9, {3.51637805487337*^9, 3.516378084716223*^9}, 3.516417861375032*^9, {3.5164179500611877`*^9, 3.516417965957616*^9}, 3.51641804834136*^9, 3.516591175169922*^9, 3.5466624175445004`*^9, {3.5466625133130684`*^9, 3.5466625344043055`*^9}, { 3.5466625914224057`*^9, 3.5466626190500546`*^9}, 3.546662845234852*^9, 3.5466628890085287`*^9, 3.5466629270569954`*^9, 3.546663960995611*^9, 3.546664111816676*^9, 3.546665023856678*^9, 3.546665308884779*^9, 3.5466653713940887`*^9, {3.5466666822911987`*^9, 3.5466666882192087`*^9}, 3.5466715630561714`*^9, 3.5466723274107137`*^9, 3.5467439335466557`*^9, 3.5467439692083187`*^9, 3.5469248487987366`*^9, 3.5469249137884502`*^9, 3.54692501004062*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ht", "=", RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{"Wt", "/", "s"}], ",", "s", ",", "t"}], "]"}]}]], "Input", CellChangeTimes->{{3.5163731554551506`*^9, 3.516373156203952*^9}, { 3.5163731948256216`*^9, 3.5163732020172343`*^9}, 3.516374013039666*^9, 3.516374301290577*^9, 3.5163744371980157`*^9, {3.516374549580613*^9, 3.51637462187114*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", "8.`"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2.`"}], " ", "t"}]], " ", RowBox[{"(", RowBox[{"0.125`", "\[VeryThinSpace]", "-", RowBox[{"0.125`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2.`", " ", "t"}]]}], "+", RowBox[{"0.25`", " ", "t"}], "+", RowBox[{"0.25`", " ", SuperscriptBox["t", "2"]}]}], ")"}]}]], "Output", CellChangeTimes->{ 3.516374493904115*^9, 3.516374631121956*^9, {3.5163747379197435`*^9, 3.516374807979467*^9}, 3.516374850161941*^9, 3.5163755002126827`*^9, 3.516375539992753*^9, 3.5163758655497246`*^9, {3.5163759490722713`*^9, 3.516375967698704*^9}, 3.5163760007551622`*^9, 3.5163765224980783`*^9, 3.516376593977404*^9, 3.516376840074237*^9, 3.5163769437136183`*^9, 3.516377037251383*^9, 3.5163772602693744`*^9, {3.516377331155899*^9, 3.516377353697939*^9}, 3.5163774943945856`*^9, 3.5163775283090453`*^9, 3.5163775886187515`*^9, 3.516377622751611*^9, {3.516377790047305*^9, 3.5163778042901297`*^9}, 3.5163779155339255`*^9, {3.5163780549045706`*^9, 3.5163780847474227`*^9}, 3.5164178616090326`*^9, {3.516417950139188*^9, 3.5164179660044155`*^9}, 3.5164180483881607`*^9, 3.516591175310547*^9, 3.546662417747301*^9, {3.546662513453469*^9, 3.546662534435506*^9}, { 3.5466625914536057`*^9, 3.5466626190812545`*^9}, 3.5466628452660513`*^9, 3.5466628890397286`*^9, 3.546662927103795*^9, 3.5466639610268116`*^9, 3.5466641118322763`*^9, 3.5466650238878784`*^9, 3.546665308915979*^9, 3.546665371425289*^9, 3.546666688250409*^9, 3.546671563087371*^9, 3.5466723274419136`*^9, 3.5467439336714563`*^9, 3.5467439692395186`*^9, 3.5469248489859366`*^9, 3.546924913928851*^9, 3.5469250100874195`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{"ht", ",", "1", ",", "0.95", ",", "1.05"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "1.1"}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.516373215277258*^9, 3.51637321657206*^9}, 3.516373247928115*^9, {3.5163736000529337`*^9, 3.5163736189289665`*^9}, { 3.516374016674472*^9, 3.5163740438965197`*^9}, {3.5163748121602745`*^9, 3.5163748452011323`*^9}, {3.54666288435972*^9, 3.546662884780921*^9}, { 3.5466679712746625`*^9, 3.546667971961064*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV1nc8lV8YAHAR2XvPe10pMiqrJOcxk0QpSlGIpCSSUalkZI+QkFFJA9Ew UrLlZyY7I8qe157J77z/3Pv5fu7zuc85z3Oe97xk22um9rQ0NDRa+IP4Tr5r y7DHLlQzXJkhKsskVqOyhWR0iuSILoCPilOknebmt2I/Eukmql9jjy/1iNVk j85/x9zpiyy2VatWK77XTC4Nl6AnBaN99zQeR+hXaspR7SI2HSKRmv5YaEpe i2Yve8q/vx3RiLwPtjvO92tKvWRSXVt/hHrODP/YXzeu6XTQ/eqyRALKCn0w sjgyr/mxtf/Fgk4SKn/wixKs+ldz7YpRz6xDKnqrxD3JokOHnm7pCV7We4aE ZaL5FAOYkG7Dw+1LHc9Rx+vxJwF3ONFo/KGyBccXaKaxjb/kFy8Ks9uwnF9P R9KZB9luzAmh3bs/rsyGv0Jbs1LvxlVKoJb1S7EzEm/QYZ/2yKu7KcijWnw3 9X0G+lXQzGE1Jo2EY1rrpnSy0MRm62IB9y5UfC7EYbLtLeK9Srfy770ispUF ugmHHMSp199meHEvol2/9X3B+x2aEZK5/ltJBbX2m4kt6L1HD88cCjqtuw95 ZbEUznV8QI5VqsLLaggZRg8zzD39iMJvarduidJCol5lJ2cdcxF3zMpjqwwd NG2V9Hxmbx4yyjh4VG1SD5XqeM5Q1/NQlrlLxztzAxQtY6pJrcxHgjmcKUPj hsiOQz5sOrwAqS0YPB9/cxRt6/6zY0qiEL17wb5g1H4cdZZ+dZ8cLUS3rfzd d9w8gTJexldMvP+M3qi1T9xQMkPG143Pj+sUIVWVOwdMZ04j0mmZt2OsX5Ee a4339/kzaO7g1vXRtq/ogl95XCKnFYpj+hw34lCC6CaUmsbrrFGup+F1N8FS FP+xWTxiiy1KGlcTmvcuRQqtvzlzdC8g5yZuhzm9MnTjYV/yngV7xJ30H+1s RzmKSdlQubhxBa2z579x0ahAm08Z6OVfXEWD99OOzTytQLfOJzr1WVxD+Q53 U6iOlShV/NPuEzdd0RklZfXp9Srk6mJhMeHsjnRekn9ftf6GpAYVTdl9PJCc EEfQVOU3tDPpvYtagif692+sbTK8GvVepHz5OXATPatNdZ2QqEF31s+7z3Tf RQvxebJ9YTWoOnHNoW+7D9J3qBtoXq1B8u33fhT/8UGTdMtmn1tqkcbiJ2m/ K75o30Fj9eDAesTmPH3nZ2EACmWxm/deqEfi/5JCNRIfoF8/b2a52DSgTunz mjX3ApG/R7r46QONaFNRPPT6mWDUlPOXdgf1O3Lfd7zkpnE44kZM59gUmhBT gg39s9MR6IB0rHqgXROaoR/J5rwUiSIW3szf+tGEFOV1yx3ropDywzZ728wf yHf4xzincAy6VytnuOd8C+rjUZj1UotHr98XSGc+akFlJVYyciPxqDlem257 fQtCdNK9IQkJSMrhdJHg/lZ0mKOm24DhCaql81eg4WlDKfnFtHmbyYj/YDf3 92/t6Hi89YPCfc8RSF2kGmy0I7+3NWlf0p8jR5bZunKlDvRJ9lwVM28aKvrJ EJD/tAOp6mYFhy2kIVuPvcvJtzrRiEC3DV1VOsrOCe52UuhC9AlvWBqfv0Hj 3ONhhjZdqKvYnTlYJAPt8DDU3Bnbha48nJ8yf5SBnmqwPBtY7UKbUwZr0mGZ KLom7OKZqm5U0C15yOvBW+QxEDl7yLIX0VGYNBzb36FUljo/Lr1+lCOb4PXk bB7SPHjb0udcP9ofq/eT+ioP9TjvUqF69qPH+stXPy3kIcHm0OH6jH7kZWO/ Ky0yHz2MNzocxPkbrXF+DqypLUB+2xvZN7t/o4ZRHY9Ms8/IAf1ImLw+gEx+ rQ7JfypBNM+XPneHDqCw+24O/dvwud8q2lP7YgC9e+t9pU2pFP1Xc1H8TfsA +qiwGANhpWjnyb/P7dUH0UH+C/c4URkavbQj69eWIVTgT+eOssrR5eg7xU1R w+j6ph+Z/LwKOQ/JDORmjyENYf0v0cN1iJnsv6Nv3wzSObWqc1OuDUVw1sk+ s55Hf/+FMWWOdiPb7pMldPVLSOSUrp/zmd/Itcrf9tqONdQyJcj53mgQiQ/q MywpraGkISbuXvNBVEfLlOENa8i67NAfWptBtB3CZ4Mt1tBtftPxA+6D6Ofn Rz4vQteQnref34XkQaT17mXKT+oaMv38zzRochBxPfmvS+fTOsrLdJWTfjCE 3ruwnhA8vIEiavokil8Oo/pEnptvzTdQdlvhNP+7YTRcKZyqbbeB1E+cCr5U OIyEhWQmnO5toEQJ5Z1LdcPIv0zPryxvAzGH5Kl/mRlG5tz3Pl6W/Icm+rvR 9L4RtPZxlufr+j+kcrq+mr5yBD27fWV89yka0ETSnT5Vo2h3KdNKoBUNmP5s sU+qH0WlW1/T912gAab550ofWkZRf/ggKdyFBhwC6hLqf48i8adWp8aCacDs 1WJH0cYoSqw0qXxWRANunt624cpjKJpNOYVbcgvQtChc25IyhnxT/x5fmNwC yuQPDs3W4+jY7k+uu+a3wO2vKZEeF8eReJnbQ9tVHK9T2cvnNI4+/xlvaqKn BY4OlhU9z3E0t/2nyVtxWlg/u+jmFD6ObLPyjl48Rgse1hUMbZ/GkVah8+HO j7SQrR77fCfrBGI3lHXk+EwLe9xVhNS4JlBP11CQfiktRDjmc2nxTyCvv5b/ 5dbTwo/GgZhDpAn0Dh059HCIFkqW/5OXUJ5ApG879AwF6ODheoPm+pkJRNPS D0W36KBpKJ7ldvoE8h8rfnbrPh1EhQ9cs8mYQEw0ybT7g+jAMmjHsnbOBOKW t6jMe0QH9Xp6d5Y/TaDtgT8O5byjA7obVo3S9RPI8EC58fNhOrCS2J7/bHYC xaSlWQaZboW5Jl8Gyf2TSPDz/a/6FltB9Zxth9vBSZTUdF6c3nor6DdsR2Va k+jlhkj//atbgUEu/rex4ST6dCr2wu2grSD9KrNN5ewk6mEJuOxcvBUKoqqj VbwnkfSNi14nZemB+iDuN++XSSSv3ncidTc9OL8UvJNaPImUaU4rjqvSQ1BK +HPp8kmkHXZ4+J4OPVT1WMbJ1kyicy/kTmZa0sPZxAxnpo5JFNc6p0gbSQ9z 2SyHts9NInqVeyM58/Rw7cRbn9HtU4h1fbl8bZUe3gzGpzDJTCHuMpcUPRoG qGFOTNohN4VIR23NulkZoGzkqM6ZvVNIw16vgkGaATx9DfwTNafQjUcsqVan GUDKNJE/2HwKDS4+Nmf5ygCLQ1X6u/2n0HMtJbfeCgZQMbEOkA+cQufDGyNz ahlAeP6m546QKdQlRV9zopMBfryJ9+WKmkI/Tlw/kDTPAByKAqHFT6ZQ8Xsj krzsNvBFhjz3Pkyhx1e3jBs/3gbpmi+Zc3qn0MnCJAZyyjawIGslGPTj9dPv o8y/2AYfihgP9P2ZQhFJzpaPP2wDGo+ICrrRKeRf193Y17ANvq7x7ZOdm0Iu MvkfXbYyAu9A+HYq/TQyHLp8J9qVEfZfNrV02jWN9O9QsnS9GEHJR1MuTX4a afP1dC3dZQTZ6um7HYrTSF3v6L6zYYyw1+12/j7laSSbrjhPecUIruurH3o1 phGL/cKlvB5GeM4tsiX36DSqH7xzsvMQEyzIXA5/fnUa/eet6hdizAScmXJn I65No0pe6nsNMyYYOCW77Ok6jYp0rTme2TLBHue657ru0yjrhXat4x0m2PZ7 kqHSexqF222D9Q9McKiAgWoZOo2MB6N2iYszw8XUJB2GVzjfrsFWSylmUExc vBzzehodcVO7+0SWGQo2JgfFM6aRAe2vJkE1/Puadsfut3h/pF0ePMeYwciP fwA+TiMVy6pSRl9mUFd33uQumUairatm80PM0O0zcyikbRolihzd3DPJDBp9 t9oH26eR0IWnr13mmCHhe+TYwc5pxD+nvz71jxkkT2unTXRNI07O2KejAizQ GDIQrtI/jeiMFCZ6D7NASECPv//4NBqvsPX57y0L1Nvo9wv+m0YSDzsvR+Sy wO3Gk1rWm9Po5Dljs5NfWEBjhEnoJQ0VFa/sl+3/jwVCWRWU5OmoKFqeq3V5 gAXarZy1ZBip6MDjkh07hFnhsK3E/i4uKgq7LPo94AErcAUu98pKUVGZWnTh kXBWUCng0dTdTkVLWxlfcMWygrkum4GlNBVZp857JT9jhS/vE84H7aQi5dZa ybwiVrA61vitSY6Keg/e9BicY4WD92OLKSpUpMDVLqZzng14JVdMenSpaDh8 2PeRPRtY0A5nfdWjomTm5ZGRK2wQpyqzLUWfili3Cn4I82KDmbjNgrMGVDS+ ZKHf/pANXJQHdtYdoaL0nl5nx0o2OLJKVbptSkWir4dKImXYQVzNR3zHeSpq lVqS+qPIDnyhGzLj2KHPGEKUVdnh0bzTiSxrKlpL3HHypzY7ZB1ZpFOwpaLO cMdRSUt2sNVFS9L2VBTjNsWZH8EOUi9NZqmXqYgZLdr0zLPD1Uu/UwbdqchL aydj/ho7+LTcQk888H50zmZHbuGAokNvVo55UlG5QdmaNgcHJLlwFxZ6UdEt 0/CYN7Ic8FBkNtLnNhVN2G+v8rDhAI/96z6DPlRUH2a2k+s7B5jcvBbBGkpF 6pFBjeNtHKDH/kHrNfbrh19uVPZwgHFvk7xOGBX5x5HLPMc5IL/vQq5HOBWJ Ny3JUWk54dKH8djWSCpq/1leZSbOCe9M77g4xVLRoakzy5JmnGBM417FnUxF Kbx+LfdOc0JHhfZ4LPbCgcycnrOcYPaf0mX+FCp6GrLuEGfLCSo+cXcFUnG9 diR1MrlwgqLnd3buZ1SUZdtbOBPCCV+9WtzHXlARZ+d57+JSTrifLnFTOouK Lm4Gnhap5ARvg6+uj7CLpN8pe1VzQiHpRQ/dWypydKeZ2tPICT/RBksfdjn3 M6v0bk54/eb+QEQOFd04+kczbIkTWOQDGNs+4P6U2205I8cFE5+7EzY/UdGM VpfwHUUueDd2rfpUIRVtKzNRfrqXC/qObw/NwVYtUXcY3scFlbc6260+U1Hs F86G63pckPPaPyzvCxWZ5BbFh57jAr8K987jxVT0LZ1399coLlDvUOHTqqCi X1Ihh/tjuGA1z2E0DHsxbdOW7jEXbDs1d74TW+r5+KPDyVxw37XwoHMlFfmm lK63veGCs6shJY+qqEjzsVP1dDn25CvV1moqyguqPEda5AINv+AcrnrcT8M6 wcAVLvjeaHLkNHYia3Pz1DoXcKevBqdg+0T16X+h5YaEB42zMg1UZBS/pmDO yQ1FaRpeGo1U9OfV7s3QXdzw3M628lgTFbFXJ6Uu23CDyjq79f5WKtoSnGZx zp4bkjcV9rpjLxhm8FRd4oadbCO577B/NhYERl/jhu59uQ3SbVT0or35mvwd bmj48/kURzs+f8OMcCGeG/h/5VCaO3A/6d37Gxu54e0wVZncTUXCEjUd6z+4 4d6n5tvHsBv3iX3f2cYNYxFfA+9hqzpVffXt5oanAd+oPdgMLXxPVEe5YeGp KeOjHjy/T/PNUrbwwIGp/MdrvVQ0eGCl9qoyD+xLG6ZN7qeieDOj8kQ1Htgi use5Ctvo2tPCanUe2Puo5uUUdu5zg9dkLR7Q6v7pc/A3ngemhIC2ozywESzW 14Ut1bEfDjrwwP3K216sA1R04frtPNZEHmjbtvnBaAivb++vPJlkHpCMTQ9z xmaeg3z9pzyQXOCmEYX9zpWhwOclDzzLNZZowf7r8vDTwgcc/+rysvkwFcVd e/Wlp44H3sh0i5waoaIap5ayrA0eMDC4ck99jIqS5FTLa2h4Qfjkt1unsF0m 48uH6XhBloFkegOb3+lcBYmZF44HeT96i217Zawylp8Xmo+HuomP43yOm9Xe irzgbTb2fhVb0WFXo5ENL8x1UcweT+Lz6Shx7a0dL/TkN+z/gB3mxMPJfokX 6i+xr9Vjj7muH//uzAsOTjzstFNUlOZd13bcmxcWU04HX8YWiL7Sa/6YF9Sn ri8qT1PRZlHmpHUDL7AXqfa+o1JRdklqeFkTL6ha0IRXY1uWxyhItvKC+OqE 2C/swurbLgNdvMByvH6AeYaK3JqPLNiP8kL2fZGtF7BHRibWL9PxAUdzciHb LBV955Zjcd/PB2y0VnTGc3g+rl4de6LBB6zaN3vOY2dXZ1eXIz5or/ue4ood 4b3Hn1OfD1Lzk/pisY2GVTcyTfkg4FzVyS7s2kIt6p8rfOAaY73dZp6KqqxP tRxL4QOhkYdC5xfwffQl/r3nMz4wsdEzdsZ+wdcVmfKCD1oO3Ll2B9u/1tJo MoMPVKHb8wm2trJtVWABHzhemnvejl267WpBcRMf/Fdz0PPwIhV9yfZ9IkfH D1HyHJ/IS/g8G5kJP2bgB+sS4xJF7ONjOxK2MPPDaQGajwexRygNce2c/HDT XNLaApsnQSDaR5wfWkT0WiOxr/hlBbXu54fd0THqq9gip9vdvV344cklmrfF y3jeFl8vDLvxQ5WLi1st9vvo227HPfmhPM1Oph37UAPJVfouP/g8VLKbwnbT vuLUFMoPbJrdAiIrVFQnR3NB6hU/vL58LcQN+/YW2WP1vfwwOF2lKLSK7ytd lk7p3/zAxajwRxL7WNDk+fuD/ND5SDlQDluKI+ea6gQ/5Mjbv0HY9aLKkU9X +IHn36Fj9tji+zQbb/AIQEycvVo2drmzqZHYYQEIMZ/pV12jIrkPSq2eRgJA 6r/IiLDjFnktm00EQH/sodQh7CveHVcCzQXg7wvtA6eweUMsQ2cvCEBLYP6A B/bF9Iu1VXcFgMtwSD4Xm7nnloFzrgAsp55h2bGO59fzDVtkgQAEe/y8Jo9t w93ZnPNZAM7/ZqtXwv5koGI1WyoAdEFlVwD7Yh7V9UaDAEg3eR62wC6LsHty a1gA/Ph4A4PXifvcePqBoCB4/hM7Noj9osf74ysRQWA1rigew/7umen1n7gg KJJOSVOxpbMZ6Zi3C0LIcef+VewW4QqB8D2C8CD55BaOv/h9Z0FNK8YQxztn bKhhD70kx6Z6C0J7X5xIAPak/+9VOh9BEDiwSzQEe8722flLftgoky8Sm0aC vGtPiCA8C3g0H48t/JhUXv5YEFi2UiyzsI8FSswMfRCEo588Rn9gFzmIGcmN CUL2FS4FwQ0qqtDrfR81KQiMcfLnRLFrKckCi1RByHE6FUTC7ugTHfi6JAh0 Bxdrd2LPnRK9ZbJVCCSqNZn3Y+80EHl9XUIIHu/k9DqNHbtTaGuhmRBInc3R j8G2PCf9auC0EMjpHQiLw5aKVTJktxQC7sNXGxKwczePRl2wFQJfsy7tp9it 7b6i7NeEQGHKfS0LmydgUvlCkBBYigQ+q8KO/l1iz/ZFCIJK9Bfmsc8INDDu KxaCXXzchUvYkke7Mm3LhMCQRe7WKvaHTwuzBdVC8GTP+dl/2M2RMndtW4VA zMsgl+kfFXFpxsQVTAlBrrdOlDh2VOLFahuSMNxJCi7Vw/b08HXzpgjD93QT VgPsc6YpEo+lhSFvg9vMEFuOud2zXk4Yfok+6zXGrrmpt1NtvzDcfapfdwqb 1kIqmNVUGA4q0Yk6YrsL/DEs8BOGTp73WcHYlvMbSz8eCMM9pSOpodi634XS JoOF4VWrWWQ4Nnfg8XVylDDodNY4PMTOXirNDEsShs9pQWvx2KNtqay2ecKw pWRPxivsM7FWjawjwqCZM3CoAlv0I7PIlXFh6I0M5a7C7vtR4FAzJQzlyWld 37DtOLi3PFgQhsiMvRdqsZ1DviltbhGBWdPd5j+wfX0UE2ZFReCk8qdvv7Df OG2xazcVgas87mPL2E6h2e+UzUXAY7X80iq2QsbZjWgLEdgwnRhaw/4wkhdn Yi0CP+90/NzALrK9XPPfVREYf+T6mm6TippOtyh8CRKB/Q6aA+zYK7rpq6nF ImCxd+kPBdv7S6pfZZkIeOmWsm3HptmbyDZWifOHrqpKY2+TiCTvrRMBu40y /53Y/Ktehys7ROB32QCbPLbSW6OE0RkRiMuQ61TBduZb2LeHIgonCyzP6GPP hU6Xm0mLwhzTL/ND2O60Y0a3ZERhaw+LqQGRn9prXaEoChr0N3QNsUNrqoPM NERhUfiVkDH26ztPOm+aiUK/E3+gGfbAkLZXeaAohIrE37yArZ3a+Fw2VBRU G6qt7bCfnj7TEB0hCvTt6fr22OfqXCUvPML5Am6wO2B3vXtaT5cmCvNjURGX sZtv/yPpFYuCbZ6NqSt2OdfnmuoFUXiRX7RwB5tcp7eguCIKrktKeXexffx/ iMevi4LotYc37mFrLo+6XaIVg0F7+hkf7M89AuJMnGIwdHHjpx/2h1fu1w13 iUGwzYPbwdjPD+4RabDB8RrHN2OIfMbMfut2YsB2eTg2FvvP+YFxmUticDEu W+YRNsX30ecHzmIgq9VzLA77xbdVC7gtBpy/ah/FY6eblMd/jBWDtuD+tSSi XjYn+BOrxaDejGZ/OvZ/1+Xu/lcrBpfZFvIJj/rRDy81iMHJWnXll9g7Xxbk nWgVg6qNSrlX2G/GRczYfovBjYwhvjfYGW6DMT5rYtC3XaImC/ttgDuXg7w4 JMxcLfyIPalh5RuwWxwEY8+Tc7HlFnTn05TEgbJSGUQ4w5a3rX+/OHBMTZnl EetDH+PP6IvD3/L6sXxiP6szEsbnxUHyu+VKIXbiVScFlYfiEMUUlldC9Evq ZOqJWHHIkjNmK8UW6jnAef2xOGzPSrAjHG/IMpedLA4v93/nKsN+tCMjb2eG OLDmVl0sx476PaIhWiEO77MvrldiB5pdOEK3KA5axv4tNUQ/THos7VfE4ah9 B6UWu+iwmXP1ujj4+eS6EZ47eOhhKK0EkF04uOuwLaV3dfBwSoBawZnD9dh7 l+dspXZJQPNYx/NG7F/xvrf0bCTA+FoVfQv2avR66Cs7CdihE2NAmDf8RjLT JeyqfyGED9+/WNrgLAGBHO/ZW7FzHQ0ZzL0lgE85mqsNO0SdO/riYwmQ7i6h 78BW6Xn2JrBBAnw+O2Z2YXMqXFrJaZKAIXbtPsKT9xQOdbZIQAJ9C3c3UV/K l8GdXRJgcY50kzDvlVaJ2mEJcMFvKj3Ya950RxhoSGCvYPRfL/ZPFtUyAyES pFTopvVjf1NOGO0RJkGH2mAj4VyrvxyuoiQo8z+0Rjgqp/xcggQJqmfOHPuN bXDC5O/YdhK8n/26SvhT4iW10L0kWFT+rDVA9EvmyduGIyQIlkl7PYS9MK+f 7neUBOmn7ZoImxTPJe03IYEK0/gyYYYThmHppiRoMdbVHyaeH96rl+9YkEB9 9mA/YdPvp3bKO5BgQiyKaRSb1Z03LcyXBM9uTCiNY19CpYna/iTY0P51knAF k1P0SgAJgjJT3QnfSqm4bxdMgv0ukfmEx6qvW2tEkeDTR1nVCWL/wj/EJpNJ kLgjQm4S+15pePyRQhKMcA4uT2GP95bGrn4mga3be+5pbPP1+ahXRSTYd+yU PGE51TPBdKUkiDsvb0u4M1P61pdvJIjmYqolvPdxydldrSTYm5gTTcUevjor zkIlgVSrItMs9vFQKZHCGRLYHfeQIFz0+pSAwxwJvo48USEcM/CVo2IR5yMF 2BAGixCa239JUDrfWkA4UZfyZ5yJDDW/o8/NYR8VMUuvlSJD/sHSR/NEPwYV Ss9Kk2G2dOsrwoNvGbsnd5BB13nPJ8LeUMTJsYsMATU6XYSz7CneJ/aQoZPx mNgCNvO7WdOeg2RwPyaZSrhaN4Jm2pwM5CTf6EVsbvZLIvdOk+HXWcZnhM91 aKlyniGDFOOtHMJLjotX9lqRQSOWXE9YOtKyw+MCGf66+W9dwvb/KZu9eY0M 4y/euRLWcq4+yxVMhlApc81l4vx2julXhpDBtT76COFcHda9nmFkyI77cpqw tNBxxt5IMugxNl8nzFrZlfs6jgw6D/leEu4UmmaDF2Ro3sfJtILtUsVb6lxC hot2MWWEZXerZZLLyPDOI6Se8ECiRVxrORk+sLl2EDZ3SXFS/0aGLXJsU4Q1 RHYKMTSQgeaHn8AqNqPrgevJXWQYvX/PkXCqqK1U/SIZSC9NaNew1ZV2Kdgu k+Hu6Vpmwq2HF9RWVsiQJLGfhzCj54MjUn/JEKwwTyHs8iPjujetJAwJd+sS hgfzpXIckvBAO8if8B+qv1X4TkkQ3qxYJ+zNYOxAkZWEc819tOvE/S0m4Fq4 SxJEXswyETY0fOM/pCAJzbeWBQh/eNGQeVBFEvo9w5UI+1nwr01pSYJ64PQl wlIVr+KMz0qCyUxfPeFGuZUXPZaSUGJe0EzYK87g4+VzktCU8aCTcJ3j2PcH NpIgZsA2SPg65y6mEgdJSAr+ska4xCr7tuINSdAOSZf+i31mJdeaM1wSxtOX vQhvtaW/lhIhCeIiv+8Qzq4zuyMXJQmhlRW+hLekLiUYxEjCsZ83wgi/1tvX 4pMgCT2sgamEF6O/6M2mS8LOALdKwlHy5bLNxbgeYW+YNoj4fWeZfEslIVdd ip1wic7CyJ5ySfgnH89NeNpCOv1hlSREsDiJED7yIETieD2u36dKOcIMfaa8 TT8l4eMJbWPCt6MGNhoWJMH8SU4Y4egn3j13liSh3MA+inDGS74v8iuSoG/D F0u4s+iQV/i6JPxZsHtCWHUsc86IlgJLYeVvCM9o3Rip46DAcNjdKsIX57f+ qJGlQEad4irhcxm6jZZyFAjIZvtL2NzGv44qTwHP4JF/hPW/033j20MBIc3I rf+Iecik/WKjRgEuzRwOwkO2NOmrOhTws1LdTtiued1LxooCv4oCjhK2DD7g UXSOAiarkyaET8JtNxNrCrxqOWZKWDd77ar7BQrM3WA/RVgqZNW2zJECFZM2 1oQHtJaNLDwocPjj1euEbd/PkYIj8f/bqcUSdvI98dfwIQV6OjziCLufyO1g jaFAddmHeMJBi+4RUXF4fZykZMLZ+1fXHydTgGPvj3TCq2Wb7S8zKfBbITqf cFQza3hlNQV0Pn5oI5x2o1HHoA7X67H1T8L5/FFrdY0UCJqm7SXcc4bnUksb Bbz8tg8S3jEgpPNngAK6pOFZwuoB3av2IxSodD2wSNhoR/K7sXEK7DN8sEL4 +hWS+OwsBSy5N/4RDmD703pjkQJ62jK0m8TzNSctdGWFAtE3j9ATzjxur+39 lwJXOi4wEi6el17d3KRA1K3rLIT/B/YxCOM= "]]}, {Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData[" 1:eJxTTMoPSmViYGAwAWIQPbcuic0wpduOAQw+2LPffqT+Vn6nPYy/8NT8wtfy J+H8des7b+fo3YLz0+0vznxT9BjOz3uq+XjLupdwPpdii/p9iw9wfp/Aaa2F CZ/h/KTbIfuZz3yD8wuPtiTlq/+C8zcW8ARLeP5FuKc6+5VBOIMDjN80/0/g lzeMcD7D5QcOe6qY4Xy1krSKEC1WOP/J1+lh3HvZ4Hyvp1m1kwo54Hy/JxO0 5eS44PxXh5MaTqzlhvP1BK/JOsfzwvlc9l8T73zmg/Pd30Z9VwoVgPNvHEph jNIRhPO3dhyJU/iK4Kexlj44d04Izk8uqt7KM0sYztdP1z7nkygC558X0uEu tRSF83eva5qtwywG51czagWcuYvgc92p8sjbIg7nP12mOGV+jQScP0VDkmVn qCScP2FW2vFEBSk4P2pK7Dme5wj+D5elP+fvk4bzHz91qjjULgPnL7I1lD6b KAvnr20tFUzXlYPz20OTvZm/IvimdxaubD8rD+fP0Jy99qy3Apxff6B3hvdO BN9XOnTpKRVFON8x73i0YCeCP18mSeXMVwRf5fDyaX7RSgj/6R7SurQPwU/7 zHLxpJYynJ+08ZNCZz+CP+EST++R4wj+vs9qP///R/ABztbzWg== "]]}, {Hue[0.1421359549995791, 0.6, 0.6], LineBox[CompressedData[" 1:eJxTTMoPSmViYGAwAWIQPbcuic0wpdsuDQze2bPffqT+Vn6nPYy/8NT8wtfy J+H8des7b+fo3YLz0+0vznxT9BjOz3uq+XjLupdwPpdii/p9iw9wfp/Aaa2F CZ/h/KTbIfuZz3yD8wuPtiTlq/+C8zcW8ARLeP5FuKc6+5VBOIMDjN80/0/g lzeMcD7D5QcOe6qY4Xy1krSKEC1WOP/J1+lh3HvZ4Hyvp1m1kwo54Hy/JxO0 5eS44PxXh5MaTqzlhvP1BK/JOsfzwvlc9l8T73zmg/Pd30Z9VwoVgPNvHEph jNIRhPO3dhyJU/iK4Kexlj44d04Izk8uqt7KM0sYztdP1z7nkygC558X0uEu tRSF83eva5qtwywG51czagWcuYvgc92p8sjbIg7nP12mOGV+jQScP0VDkmVn qCScP2FW2vFEBSk4P2pK7Dme5wj+D5elP+fvk4bzHz91qjjULgPnL7I1lD6b KAvnr20tFUzXlYPz20OTvZm/IvimdxaubD8rD+fP0Jy99qy3Apxff6B3hvdO BN9XOnTpKRVFON8x73i0YCeCP18mSeXMVwRf5fDyaX7RSgj/6R7SurQPwU/7 zHLxpJYynJ+08ZNCZz+CP+EST++R4wj+vs9qP///R/AB/zZvUQ== "]]}, {Hue[0.37820393249936934`, 0.6, 0.6], LineBox[CompressedData[" 1:eJxTTMoPSmViYGAwAWIQPbcuic0wpdvu7BkQ+GDPfvuR+lv5nfYw/sJT8wtf y5+E89et77ydo3cLzk+3vzjzTdFjOD/vqebjLetewvlcii3q9y0+wPl9Aqe1 FiZ8hvOTbofsZz7zDc4vPNqSlK/+C87fWMATLOH5F+Ge6uxXBuEMDjB+0/w/ gV/eMML5DJcfOOypYobz1UrSKkK0WOH8J1+nh3HvZYPzvZ5m1U4q5IDz/Z5M 0JaT44LzXx1OajixlhvO1xO8Jusczwvnc9l/TbzzmQ/Od38b9V0pVADOv3Eo hTFKRxDO39pxJE7hK4Kfxlr64Nw5ITg/uah6K88sYThfP137nE+iCJx/XkiH u9RSFM7fva5ptg6zGJxfzagVcOYugs91p8ojb4s4nP90meKU+TUScP4UDUmW naGScP6EWWnHExWk4PyoKbHneJ4j+D9clv6cv08azn/81KniULsMnL/I1lD6 bKIsnL+2tVQwXVcOzm8PTfZm/orgm95ZuLL9rDycP0Nz9tqz3gpwfv2B3hne OxF8X+nQpadUFOF8x7zj0YKdCP58mSSVM18RfJXDy6f5RSsh/Kd7SOvSPgQ/ 7TPLxZNaynB+0sZPCp39CP6ESzy9R44j+Ps+q/38/x/BBwDRVOw9 "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, ImageSize->{533., Automatic}, PlotRange->{{0, 10}, {0, 1.1}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Automatic}]], "Output", CellChangeTimes->{ 3.516374493935315*^9, 3.5163746311531563`*^9, {3.516374737950944*^9, 3.516374808026267*^9}, 3.5163748502243414`*^9, 3.516375500275083*^9, 3.516375540055153*^9, 3.5163758656121244`*^9, {3.516375949119071*^9, 3.516375967745504*^9}, 3.516376000801962*^9, 3.516376522560478*^9, 3.5163765940242043`*^9, 3.5163768401366367`*^9, 3.516376943776019*^9, 3.5163770373137827`*^9, 3.5163772603161745`*^9, {3.516377331202699*^9, 3.5163773537447386`*^9}, 3.5163774944413857`*^9, 3.5163775283558455`*^9, 3.516377588665551*^9, 3.5163776227984114`*^9, {3.5163777900941052`*^9, 3.5163778043369303`*^9}, 3.516377915596326*^9, {3.5163780549669704`*^9, 3.516378084809823*^9}, 3.5164178616714325`*^9, {3.516417950185988*^9, 3.516417966051216*^9}, 3.5164180484505606`*^9, 3.516591175419922*^9, 3.5466624177941008`*^9, {3.5466625134846687`*^9, 3.5466625344823055`*^9}, { 3.546662591484806*^9, 3.5466626191124544`*^9}, 3.5466628452972517`*^9, 3.5466628890709286`*^9, 3.5466629271349955`*^9, 3.5466639610580115`*^9, 3.5466641118634768`*^9, 3.5466650239190784`*^9, 3.546665308947179*^9, 3.5466653714564886`*^9, 3.5466666882816095`*^9, 3.5466679729126654`*^9, 3.5466715631185713`*^9, 3.546672327473114*^9, 3.5467439337026563`*^9, 3.5467439692707186`*^9, 3.5469248490639367`*^9, 3.5469249140068507`*^9, 3.5469250101342196`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"hf", "[", "t_", "]"}], "=", RowBox[{"Abs", "[", "ht", "]"}]}]], "Input", CellChangeTimes->{{3.516374890534812*^9, 3.5163749181312604`*^9}}], Cell[BoxData[ RowBox[{"8.`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2.`"}], " ", RowBox[{"Re", "[", "t", "]"}]}]], " ", RowBox[{"Abs", "[", RowBox[{"0.125`", "\[VeryThinSpace]", "-", RowBox[{"0.125`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2.`", " ", "t"}]]}], "+", RowBox[{"0.25`", " ", "t"}], "+", RowBox[{"0.25`", " ", SuperscriptBox["t", "2"]}]}], "]"}]}]], "Output", CellChangeTimes->{ 3.516375500321883*^9, 3.5163755401019526`*^9, 3.516375865658925*^9, { 3.5163759491658716`*^9, 3.516375967792304*^9}, 3.5163760008487625`*^9, 3.516376522607279*^9, 3.516376594071004*^9, 3.516376840183437*^9, 3.5163769438228188`*^9, 3.5163770373605833`*^9, 3.516377260362975*^9, { 3.5163773312494993`*^9, 3.5163773537915387`*^9}, 3.516377494488186*^9, 3.5163775284026456`*^9, 3.5163775887123513`*^9, 3.5163776228452115`*^9, { 3.5163777901409054`*^9, 3.5163778043837304`*^9}, 3.5163779156431255`*^9, { 3.51637805501377*^9, 3.5163780848566227`*^9}, 3.516417861702633*^9, { 3.516417950232788*^9, 3.516417966098016*^9}, 3.516418048481761*^9, 3.516591175482422*^9, 3.5466624178253007`*^9, {3.546662513500269*^9, 3.5466625344979057`*^9}, {3.5466625915004063`*^9, 3.5466626191124544`*^9}, 3.5466628452972517`*^9, 3.546662889102129*^9, 3.5466629271505957`*^9, 3.5466639610736113`*^9, 3.546664111894676*^9, 3.546665023934678*^9, 3.546665308962779*^9, 3.546665371472089*^9, 3.5466666882972093`*^9, 3.5466715631497717`*^9, 3.546672327488714*^9, 3.546743933718256*^9, 3.5467439692707186`*^9, 3.546924849110737*^9, 3.546924914053651*^9, 3.54692501018102*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:0412\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:0438\ \:044f", " ", "\:0438", " ", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\ \:0432\:0430\:043d\:0438\:0435"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"k", "=", "10000"}], ";"}], RowBox[{"(*", RowBox[{ RowBox[{"\:043a\:043e\:043b", "-", RowBox[{ "\:0432\:043e", " ", "\:043f\:0440\:043e\:043c\:0435\:0436\:0443\:0442\:043a\:043e\:0432"}]}\ ], ",", " ", RowBox[{ "\:0433\:0434\:0435", " ", "\:043d\:0430\:0445\:043e\:0434\:0438\:043c", " ", "\:0440\:0435\:0448\:0435\:043d\:0438\:0435"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dt", "=", RowBox[{"10", "/", "k"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hust", "=", RowBox[{"hf", "[", "10000", "]"}]}], RowBox[{"(*", RowBox[{ "\:0443\:0441\:0442\:0430\:043d\:043e\:0432\:0438\:0432\:0448\:0435\:0435\ \:0441\:044f", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hmax", "=", RowBox[{"hf", "[", "0", "]"}]}], RowBox[{"(*", RowBox[{ "\:0437\:0430\:0434\:0430\:0434\:0438\:043c", " ", "\:043d\:0430\:0447\:0430\:043b\:044c\:043d\:043e\:0435", " ", RowBox[{ "\:043c\:0430\:043a\:0441", ".", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}]}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "\:0427\:0438\:0441\:043b\:0435\:043d\:043d\:043e\:0435", " ", "\:0440\:0435\:0448\:0435\:043d\:0438\:0435", " ", "\:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:0439", " ", "h", RowBox[{"(", "t", ")"}]}], "=", RowBox[{ RowBox[{"0.95", " ", "\:0438", " ", "h", RowBox[{"(", "t", ")"}]}], "=", "1.05"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{ "\:043e\:0431\:044a\:0435\:0434\:0438\:043d\:0435\:043d\:043d\:044b\ \:0445", " ", "\:0432", " ", "\:043e\:0434\:043d\:043e"}], " ", "|", RowBox[{ RowBox[{"h", RowBox[{"(", "t", ")"}]}], "-", "1"}], "|", RowBox[{"-", "0.05"}]}], "=", RowBox[{ "0", " ", "\:0438", " ", "\:043e\:043f\:0440\:0435\:0434\:0435\:043b\:0435\:043d\:0438\:0435", " ", "\:043c\:0430\:043a\:0441\:0438\:043c\:0443\:043c\:0430"}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"k", "+", "1"}]}], ",", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"hf", "[", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}]}], "]"}], "-", "hust"}], "]"}], "-", "0.05"}], "]"}], "\[NotEqual]", RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"hf", "[", RowBox[{"dt", "*", RowBox[{"(", "i", ")"}]}], "]"}], "-", "hust"}], "]"}], "-", "0.05"}], "]"}]}], ",", RowBox[{"Treg", "=", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"i", "-", "0.5"}], ")"}]}]}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"hmax", "<", RowBox[{"hf", "[", RowBox[{"dt", "*", "i"}], "]"}]}], ",", RowBox[{"hmax", "=", RowBox[{"hf", "[", RowBox[{"dt", "*", "i"}], "]"}]}]}], "]"}]}], "}"}], ";", RowBox[{"i", "++"}]}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sigmaper", "=", RowBox[{"Round", "[", RowBox[{ RowBox[{"100", "*", RowBox[{ RowBox[{"(", RowBox[{"hmax", "-", "hust"}], ")"}], "/", "hust"}]}], ",", "0.01"}], "]"}]}], ";"}], " ", RowBox[{ "(*", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\:0440\ \:043e\:0432\:0430\:043d\:0438\:0435", "*)"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "Treg", ",", "\"\< \:0441\>\""}], "]"}], " ", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "sigmaper", ",", "\"\< %\>\""}], "]"}], " "}]}]], "Input", CellChangeTimes->{{3.5163749539645233`*^9, 3.516374954245324*^9}, { 3.5163749848681774`*^9, 3.51637529920793*^9}, {3.5163753388475995`*^9, 3.516375545998763*^9}, {3.5164180254405203`*^9, 3.516418028732126*^9}, { 3.5466653012719655`*^9, 3.5466653033467693`*^9}, {3.5466653634848747`*^9, 3.5466653674784813`*^9}, {3.5466724043656487`*^9, 3.5466724071112537`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Treg = \"\>", "\[InvisibleSpace]", "3.1475`", "\[InvisibleSpace]", "\<\" \:0441\"\>"}], SequenceForm["Treg = ", 3.1475, " \:0441"], Editable->False]], "Print", CellChangeTimes->{ 3.5163755012266846`*^9, 3.5163755409911547`*^9, 3.516375866548126*^9, { 3.516375950055073*^9, 3.5163759686815057`*^9}, 3.516376001753564*^9, 3.51637652349648*^9, 3.5163765949758053`*^9, 3.5163768410726385`*^9, 3.5163769447120204`*^9, 3.5163770382653847`*^9, 3.5163772612521763`*^9, { 3.5163773321387005`*^9, 3.5163773546963406`*^9}, 3.5163774953773875`*^9, 3.5163775293074474`*^9, 3.516377589617153*^9, 3.5163776237500134`*^9, { 3.5163777910301065`*^9, 3.5163778052885323`*^9}, 3.516377916532327*^9, { 3.516378055902972*^9, 3.5163780857458243`*^9}, 3.5164178627010345`*^9, { 3.51641795116879*^9, 3.516417967049618*^9}, 3.5164180493865623`*^9, 3.516591176216797*^9, 3.546662418324502*^9, {3.54666251409307*^9, 3.5466625350439067`*^9}, {3.546662592077607*^9, 3.5466626196272554`*^9}, 3.5466628458744526`*^9, 3.5466628896637297`*^9, 3.5466629277277966`*^9, 3.546663961572812*^9, 3.5466641124250774`*^9, 3.546665024449479*^9, 3.54666530947758*^9, 3.54666537198689*^9, 3.54666668881201*^9, 3.5466715636645727`*^9, 3.546672328003515*^9, 3.546743934248657*^9, 3.5467439697855196`*^9, 3.5469248497347383`*^9, 3.546924915457653*^9, 3.5469250111170216`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Sigma = \"\>", "\[InvisibleSpace]", "0.`", "\[InvisibleSpace]", "\<\" %\"\>"}], SequenceForm["Sigma = ", 0., " %"], Editable->False]], "Print", CellChangeTimes->{ 3.5163755012266846`*^9, 3.5163755409911547`*^9, 3.516375866548126*^9, { 3.516375950055073*^9, 3.5163759686815057`*^9}, 3.516376001753564*^9, 3.51637652349648*^9, 3.5163765949758053`*^9, 3.5163768410726385`*^9, 3.5163769447120204`*^9, 3.5163770382653847`*^9, 3.5163772612521763`*^9, { 3.5163773321387005`*^9, 3.5163773546963406`*^9}, 3.5163774953773875`*^9, 3.5163775293074474`*^9, 3.516377589617153*^9, 3.5163776237500134`*^9, { 3.5163777910301065`*^9, 3.5163778052885323`*^9}, 3.516377916532327*^9, { 3.516378055902972*^9, 3.5163780857458243`*^9}, 3.5164178627010345`*^9, { 3.51641795116879*^9, 3.516417967049618*^9}, 3.5164180493865623`*^9, 3.516591176216797*^9, 3.546662418324502*^9, {3.54666251409307*^9, 3.5466625350439067`*^9}, {3.546662592077607*^9, 3.5466626196272554`*^9}, 3.5466628458744526`*^9, 3.5466628896637297`*^9, 3.5466629277277966`*^9, 3.546663961572812*^9, 3.5466641124250774`*^9, 3.546665024449479*^9, 3.54666530947758*^9, 3.54666537198689*^9, 3.54666668881201*^9, 3.5466715636645727`*^9, 3.546672328003515*^9, 3.546743934248657*^9, 3.5467439697855196`*^9, 3.5469248497347383`*^9, 3.546924915457653*^9, 3.5469250111170216`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"\:041a\:043e\:043b", "-", RowBox[{ "\:0432\:043e", " ", "\:043f\:0435\:0440\:0438\:043e\:0434\:043e\:0432", " ", "\:043a\:0432\:0430\:043d\:0442\:043e\:0432\:0430\:043d\:0438\:044f\ ", " ", "\:0437\:0430", " ", "\:0432\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:0438\ \:044f"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"Nkvant", "=", RowBox[{"Ceiling", "[", RowBox[{"Treg", "/", "Tkvant"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.5163756536545525`*^9, 3.5163756544657536`*^9}, { 3.5163757406559052`*^9, 3.5163758443336873`*^9}, {3.516375944875864*^9, 3.516375945031864*^9}, {3.5163759898039427`*^9, 3.5163759960907536`*^9}}], Cell[BoxData["3"], "Output", CellChangeTimes->{ 3.5163758667041264`*^9, {3.5163759502422733`*^9, 3.5163759688687057`*^9}, 3.516376001893964*^9, 3.5163765236680803`*^9, 3.516376595163006*^9, 3.5163768412286386`*^9, 3.5163769448524203`*^9, 3.5163770384057846`*^9, 3.516377261392576*^9, {3.516377332279101*^9, 3.5163773548523407`*^9}, 3.516377495517788*^9, 3.516377529432247*^9, 3.516377589741953*^9, 3.5163776238904133`*^9, {3.5163777911705074`*^9, 3.5163778053665323`*^9}, 3.5163779166415277`*^9, {3.516378056043372*^9, 3.5163780858706245`*^9}, 3.516417862747834*^9, {3.5164179512779903`*^9, 3.516417967158818*^9}, 3.5164180494957623`*^9, 3.516591176294922*^9, 3.546662418355702*^9, { 3.54666251410867*^9, 3.546662535059507*^9}, {3.546662592108807*^9, 3.5466626196428556`*^9}, 3.546662845905653*^9, 3.5466628896949296`*^9, 3.5466629277589965`*^9, 3.5466639616040125`*^9, 3.5466641124406776`*^9, 3.5466650244650793`*^9, 3.5466653094931803`*^9, 3.5466653720024896`*^9, 3.54666668884321*^9, 3.5466715636801724`*^9, 3.5466723280191145`*^9, 3.546743934279857*^9, 3.5467439698167195`*^9, 3.546924849781538*^9, 3.5469249155044537`*^9, 3.5469250111482215`*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"Stextrop", "=", "3"}], RowBox[{"(*", RowBox[{ "\:0421\:0442\:0435\:043f\:0435\:043d\:044c", " ", "\:043f\:043e\:043b\:0438\:043d\:043e\:043c\:0430", " ", "\:044d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:044f\:0446\:0438\ \:0438"}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Masotklon", "=", RowBox[{"{", "}"}]}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0441\:0441\:0438\:0432", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:0439"}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", "Stextrop", "}"}]}], "]"}]}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0441\:0441\:0438\:0432", " ", "\:043f\:043e\:0441\:043b\:0435\:0434\:043d\:0438\:0445", " ", "\:0438\:0437\:043c\:0435\:0440\:0435\:043d\:043d\:044b\:0445", " ", RowBox[{"(", "\:0440\:0435\:0430\:043b\:044c\:043d\:044b\:0445", ")"}], " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0439", " ", "\:0438\:0441\:043f\:043e\:043b\:044c\:0437\:0443\:0435\:043c\:044b\:0445\ ", " ", "\:043f\:0440\:0438", " ", "\:044d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:044f\:0446\:0438\ \:0438"}], "*)"}], ";"}]}], "Input", CellChangeTimes->{{3.5163760323296175`*^9, 3.5163761120613575`*^9}, { 3.5163761944919024`*^9, 3.51637629556448*^9}, {3.516376478911602*^9, 3.516376482858409*^9}, {3.5163771007434945`*^9, 3.5163771481363773`*^9}, 3.516377436658884*^9}], Cell[BoxData[ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", RowBox[{"Stextrop", "-", "1"}]}], ",", RowBox[{"i", "<=", "Nkvant"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Ttek", "=", RowBox[{"Tkvant", "*", "i"}]}], ";", "\[IndentingNewLine]", RowBox[{"hTtek", "=", RowBox[{"hf", "[", "Ttek", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "1", "]"}], "]"}], "=", RowBox[{"hf", "[", RowBox[{"Ttek", "-", RowBox[{"Tkvant", "*", "2"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "2", "]"}], "]"}], "=", RowBox[{"hf", "[", RowBox[{"Ttek", "-", RowBox[{"Tkvant", "*", "1"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "3", "]"}], "]"}], "=", RowBox[{"hf", "[", RowBox[{"Ttek", "-", RowBox[{"Tkvant", "*", "0"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "\:042d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:044f\:0446\ \:0438\:044f", " ", "\:041c\:0435\:0442\:043e\:0434", " ", "\:0410\:0434\:0430\:043c\:0441\:0430"}], "*)"}], "\[IndentingNewLine]", RowBox[{"hextrop", "=", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "3", "]"}], "]"}], "+", RowBox[{"Tkvant", "*", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"23", "*", RowBox[{"hizm", "[", RowBox[{"[", "3", "]"}], "]"}]}], "-", RowBox[{"16", "*", RowBox[{"hizm", "[", RowBox[{"[", "2", "]"}], "]"}]}], "+", RowBox[{"5", "*", RowBox[{"hizm", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ")"}], "/", "12"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Masotklon", "=", RowBox[{"Append", "[", RowBox[{"Masotklon", ",", RowBox[{"Abs", "[", RowBox[{ RowBox[{"hf", "[", RowBox[{"Ttek", "+", "Tkvant"}], "]"}], "-", "hextrop"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "}"}]}], "]"}], ";"}]], "Input", CellChangeTimes->{{3.5163763045188956`*^9, 3.5163763836110344`*^9}, { 3.5163764362767267`*^9, 3.5163764572119637`*^9}, {3.516376577332175*^9, 3.5163765815909824`*^9}, {3.516376919034375*^9, 3.5163769345252028`*^9}, { 3.516377230988123*^9, 3.5163772464633503`*^9}, 3.5163773480819287`*^9, { 3.5163774539125147`*^9, 3.5163774570481205`*^9}, {3.516377487109373*^9, 3.516377487624174*^9}, 3.5163775230050364`*^9, {3.5163775575902967`*^9, 3.516377573861125*^9}, 3.516377612876794*^9, {3.516377799344921*^9, 3.5163777997973223`*^9}, {3.516377893756287*^9, 3.516377909403115*^9}, { 3.5163779812568407`*^9, 3.516378015920102*^9}, {3.516378050988963*^9, 3.516378064888588*^9}, {3.5164179455059795`*^9, 3.516417961870408*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0435\:043c\:0430\:0442\:0438\:0447\:0435\:0441\:043a\ \:043e\:0435", " ", "\:043e\:0436\:0438\:0434\:0430\:043d\:0438\:0435", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:044f", " ", "\:0438\:0441\:0442\:0438\:043d\:043d\:043e\:0433\:043e", " ", "\:0443\:0433\:043b\:0430", " ", "\:043a\:0440\:0435\:043d\:0430", " ", "\:043e\:0442", " ", "\:044d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:0438\:0440\:043e\ \:0432\:0430\:043d\:043d\:043e\:0433\:043e", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Matozh", "=", RowBox[{"Mean", "[", "Masotklon", "]"}], " "}]}]], "Input", CellChangeTimes->{{3.5163766231806555`*^9, 3.516376754570486*^9}, { 3.5163770222129564`*^9, 3.5163770244125605`*^9}}], Cell[BoxData["1.2897357104405227`"], "Output", CellChangeTimes->{ 3.5163768412910385`*^9, 3.5163769449148207`*^9, 3.516377038468185*^9, 3.5163772614549766`*^9, {3.516377332341501*^9, 3.5163773549147406`*^9}, 3.5163774958141885`*^9, 3.5163775296038475`*^9, 3.516377589804353*^9, 3.516377623952813*^9, {3.5163777912329073`*^9, 3.516377805428932*^9}, 3.5163779167039275`*^9, {3.5163780561057725`*^9, 3.516378085964225*^9}, 3.5164178628414345`*^9, {3.51641795134039*^9, 3.516417967221218*^9}, 3.5164180495893626`*^9, 3.516591176357422*^9, 3.5466624184181023`*^9, { 3.5466625141398697`*^9, 3.546662535090707*^9}, {3.546662592140007*^9, 3.5466626196740556`*^9}, 3.5466628459212527`*^9, 3.5466628897105293`*^9, 3.5466629277745967`*^9, 3.5466639616196127`*^9, 3.546664112471877*^9, 3.5466650244962792`*^9, 3.5466653095243797`*^9, 3.5466653720336895`*^9, 3.54666668887441*^9, 3.5466715637113724`*^9, 3.546672328050315*^9, 3.546743934311057*^9, 3.5467439698323193`*^9, 3.5469248498751383`*^9, 3.5469249155512533`*^9, 3.5469250111950216`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "\:0421\:0440\:0435\:0434\:043d\:0435\:0435", " ", "\:043a\:0432\:0430\:0434\:0440\:0430\:0442\:0438\:0447\:0435\:0441\:043a\ \:043e\:0435", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:0435", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{"sigma", "=", RowBox[{"StandardDeviation", "[", "Masotklon", "]"}], " "}]}]], "Input", CellChangeTimes->{{3.51637672255923*^9, 3.516376738954859*^9}, 3.516377323402685*^9}], Cell[BoxData["0.10780736503093923`"], "Output", CellChangeTimes->{ 3.516376841322239*^9, 3.516376944946021*^9, 3.516377038499385*^9, 3.5163772614861765`*^9, {3.5163773324975014`*^9, 3.516377354945941*^9}, 3.516377495970188*^9, 3.5163775296662474`*^9, 3.5163775898511534`*^9, 3.5163776239840136`*^9, {3.516377791264107*^9, 3.516377805475732*^9}, 3.5163779167351274`*^9, {3.5163780561369724`*^9, 3.5163780860110245`*^9}, 3.5164178633718357`*^9, {3.51641795138719*^9, 3.516417967268018*^9}, 3.5164180496205626`*^9, 3.516591176419922*^9, 3.546662418558502*^9, { 3.54666251415547*^9, 3.546662535106307*^9}, {3.5466625921712074`*^9, 3.5466626196896553`*^9}, 3.5466628459524527`*^9, 3.5466628897417297`*^9, 3.5466629277901964`*^9, 3.546663961650812*^9, 3.5466641124874773`*^9, 3.546665024511879*^9, 3.54666530953998*^9, 3.5466653720492897`*^9, 3.54666668889001*^9, 3.5466715637269726`*^9, 3.546672328065915*^9, 3.5467439343266573`*^9, 3.5467439698479195`*^9, 3.546924850218339*^9, 3.5469249155824537`*^9, 3.546925011210622*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:041e\:0420\:041e\:0413", " ", "\:0421\:0420\:0410\:0411\:0410\:0422\:042b\:0412\:0410\:041d\:0418\:042f\ "}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"b", "=", "3"}], " ", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"1", ".."}], "5"}], " ", "-", " ", RowBox[{ "\:0434\:043e\:0432\:0435\:0440\:0438\:0442\:0435\:043b\:044c\:043d\ \:0430\:044f", " ", "\:0433\:0440\:0430\:043d\:0438\:0446\:0430"}]}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{"Porog", "=", RowBox[{"Matozh", "+", RowBox[{"b", "*", "sigma", " "}]}]}]}]}]], "Input", CellChangeTimes->{{3.51637675677009*^9, 3.5163768082813807`*^9}, { 3.516376860369872*^9, 3.5163768651454806`*^9}, {3.5163770298881702`*^9, 3.516377030683771*^9}, 3.516377764806461*^9}], Cell[BoxData["1.6131578055333404`"], "Output", CellChangeTimes->{ 3.5163768413690386`*^9, 3.5163769449928207`*^9, 3.5163770385305853`*^9, 3.516377261532977*^9, {3.5163773325443015`*^9, 3.516377354977141*^9}, 3.5163774960169888`*^9, 3.516377529697448*^9, 3.516377589882354*^9, 3.5163776239996133`*^9, {3.516377791295307*^9, 3.5163778055069323`*^9}, 3.516377916766328*^9, {3.516378056183772*^9, 3.516378086042225*^9}, 3.5164178634030356`*^9, {3.5164179514183903`*^9, 3.516417967299218*^9}, 3.5164180496361628`*^9, 3.516591176451172*^9, 3.546662418620902*^9, { 3.54666251417107*^9, 3.5466625351219068`*^9}, {3.546662592186807*^9, 3.5466626197052555`*^9}, 3.546662845968053*^9, 3.54666288975733*^9, 3.5466629278057966`*^9, 3.5466639616664124`*^9, 3.5466641125030775`*^9, 3.546665024527479*^9, 3.54666530955558*^9, 3.54666537206489*^9, 3.54666668890561*^9, 3.5466715637425723`*^9, 3.546672328081515*^9, 3.546743934357857*^9, 3.5467439698635197`*^9, 3.546924850358739*^9, 3.5469249156136537`*^9, 3.5469250112262216`*^9}] }, Open ]] }, WindowSize->{1078, 738}, WindowMargins->{{11, Automatic}, {Automatic, 7}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (November 7, 2010)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 2698, 81, 292, "Input"], Cell[3246, 103, 521, 13, 52, "Input"], Cell[CellGroupData[{ Cell[3792, 120, 1744, 46, 112, "Input"], Cell[5539, 168, 1401, 20, 30, "Output"], Cell[6943, 190, 1383, 19, 30, "Output"], Cell[8329, 211, 1371, 19, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9737, 235, 2280, 59, 92, "Input"], Cell[12020, 296, 2139, 42, 81, "Output"], Cell[14162, 340, 2206, 43, 81, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16405, 388, 677, 15, 52, "Input"], Cell[17085, 405, 1550, 26, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18672, 436, 239, 6, 31, "Input"], Cell[18914, 444, 1618, 28, 57, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20569, 477, 573, 13, 52, "Input"], Cell[21145, 492, 1684, 30, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[22866, 527, 252, 6, 31, "Input"], Cell[23121, 535, 1657, 29, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[24815, 569, 700, 14, 52, "Input"], Cell[25518, 585, 4658, 79, 371, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30213, 669, 352, 8, 31, "Input"], Cell[30568, 679, 114, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30719, 685, 227, 5, 31, InheritFromParent], Cell[30949, 692, 109, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31095, 698, 1322, 33, 92, "Input"], Cell[32420, 733, 285, 4, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32742, 742, 382, 8, 31, "Input"], Cell[33127, 752, 1392, 19, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34556, 776, 515, 11, 52, "Input"], Cell[35074, 789, 1413, 20, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[36524, 814, 403, 8, 31, "Input"], Cell[36930, 824, 1796, 32, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38763, 861, 661, 14, 31, "Input"], Cell[39427, 877, 15353, 257, 356, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[54817, 1139, 178, 4, 31, "Input"], Cell[54998, 1145, 1686, 31, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[56721, 1181, 5099, 133, 232, "Input"], Cell[CellGroupData[{ Cell[61845, 1318, 1445, 23, 23, "Print"], Cell[63293, 1343, 1429, 23, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[64771, 1372, 777, 15, 52, "Input"], Cell[65551, 1389, 1213, 17, 30, "Output"] }, Open ]], Cell[66779, 1409, 1563, 36, 72, "Input"], Cell[68345, 1447, 3120, 74, 212, "Input"], Cell[CellGroupData[{ Cell[71490, 1525, 886, 17, 52, "Input"], Cell[72379, 1544, 1085, 15, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[73501, 1564, 578, 13, 52, "Input"], Cell[74082, 1579, 1074, 15, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[75193, 1599, 886, 23, 72, "Input"], Cell[76082, 1624, 1068, 15, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)